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THE  MECHANICAL  PROPERTIES  OF  SEA  ICE 

by 

W.F.  Weeks  and  A.  Assur 
I.  INTRODUCTION 

Although  sea  ice  covers  up  to  12%  of  the  surface  of  the  oceans,  it  is 
commonly  considered  somewhat  of  a  scientific  oddity.  It  is  still  rare,  ex¬ 
cluding  glaciological  circles,  to  find  a  person  who  has  even  seen  sea  ice. 
Therefore,  it  is  not  very  surprising  that  careful,  detailed  observations  such 
as  those  of  Malmgren,7*r  concerning  the  physical  characteristics  of  this  ma¬ 
terial  have  been  few.  The  most  notable  exceptions  to  this  are  studies  of  the 
mechanical  properties  of  sea  ice.  The  need  for  reliable  estimates  of  the 
strength  of  sea  ice  was  clear  to  Russian  investigators  in  the  lo90,s.*®,,,  The 
motivation  for  this  interest  was  ice-breaker  design  and  later  over -ice  trans¬ 
port  related  to  the  effective  operation  of  the  Northern  Sea  Route.  Finnish 
and  Japanese  investigators  became  interested  in  the  late  1940's  because  of 
pack  ice  problems  in  the  Baltic  and  Okhotsk  Seas.  The  Canadians  and 
Americans  became  involved  in  the  mid-fifties  as  the  result  of  problems 
associated  with  ship  and  air  resupply  in  the  Arctic.  Because  of  linguistic 
barriers  and  the  fact  that  the  publication  of  research  results  has  been  both 
sporadic  and  scattered  through  a  wide  array  of  different  types  of  technical 
journals,  the  study  of  the  properties  of  sea  ice  has  proceeded  in  a  relatively 
independent  manner  in  these  different  countries.  It  is  hoped  that  this  review 
will  provide  the  reader  with  an  estimate  of  the  state  of  thinking  of  each  of 
the  main  groups  investigating  sea  ice  and  an  overall  appraisal  of  the  field  as 
a  whole. 

In  the  last  10  years,  because  of  the  pronounced  increase  in  activities  in 
the  Arctic  and  Antarctic  there  has  been  a  steady  demand  for  the  solution  of 
increasingly  difficult  problems  that  depend  upon  a  knowledge  of  the  mechan¬ 
ical  properties  of  sea  ice:  the  design  of  a  barge  that  will  withstand  the 
forces  exerted  by  heavy  arctic  pack  ice,  the  design  of  a  new  generation  of 
ice  breakers,  the  design  of  off-shore  drilling  platforms  that  will  not  be  de¬ 
stroyed  by  rapidly  moving  ice  floes,  the  prediction  of  the  effects  of  the  long 
term  loading  of  an  ice  sheet,  the  calculation  of  the  forces  exerted  on  piers 
as  an  ice  sheet  is  raised  by  the  tide.  Other  problems  are  traffic  over  sea 
ice,  unloading  on  the  edge  of  an  ice  sheet,  and  the  use  of  ice  runways  for 
logistical  operations. 

Fortunately,  most  of  these  problems  have  certain  features  in  common. 
For  example,  consider  the  problem  of  calculating  the  failure  of  an  ice  sheet 
under  either  bending  or  buckling 


Mz 


<r  -  - 


(1.1) 


where  a  is  the  tensile  stress  exerted  on  the  extreme  ice  fiber,  M  is  the 
bending  moment,  I  is  the  moment  of  inertia  about  the  neutral  axis  and  ze 
is  the  distance  from  the  neutral  ads  to  the  extreme  fiber  of  the  ice.  In 
Cartesian  coordinates 


♦Superscript  numbers  refer  to  cited  literature  (p.  67  to  80). 
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and  in  polar  coordinates  with  symmetry 


where  the  deflection  w  is  given  by 


w 


Pm 

kf2 


(1.2) 


(1.3) 


(1.4) 


Here  xj  are  the  coordinates  of  the  system  (i  =  1,  2),  p  is  Poisson's  ratio, 
k  is  the  foundation  modulus  (in  this  case  the. density  of  sea  water),  and  P  is 
the  applied  force.  The  action  radius  f  is  given  by 


I  = 


(1.5) 


where  D  is  the  flexural  rigidity 


D 


ft  E(z)  •  (z  -  zR)2 

J  - 1 - 
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(1.6) 


z  is  the  vertical  distance  from  the  top,  zn  is  the  distance  from  the  top  to  the 
neutral  axis,  and  E(z)  is  the  elastic  or  Young's  modulus  which  varies  through 
out  the  ice  thickness.  A  method  for  calculating  these  values  is  given  at  the 
end  of  this  paper.  The  value  of  u  is  determined  by 


w  = 


a. 


,  f  l  bi  Ci 

£  Vf  ’  T*  T’ 


(1.7) 


where  a^,  bj  and  cj  are  lengths  associated  with  where  the  force  is  applied, 
the  distribution  of  the  force,  and  the  boundary  conditions  respectively. 
Therefore  at  failure 


(1.8) 


where  o-£  is  the  failure  strength  of  the  ice. 

For  the  case  of  a  homogeneous,  isotropic  ice  sheet,  the  neutral  axis 
is  in  the  center  of  the  sheet  and 


(1.9) 
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Also  in  this  case 
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and 
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(1.11) 


The  force  required  for  failure  is,  therefore, 

-  ^ 

13  6  f(u ,  (jl) 


(1.12) 


which  in  the  case  of  a  circular  load  with  radius  a  acting  on  an  infinite  sheet 
gives195 


(1.13) 


where  kci'  indicates  a  Kelvin  function. 

Therefore,  the  material  parameters  or  their  viscous  analogs  that  are 
needed  for  the  practical  solution  of  this  type  of  problem  are  E,  p  and  a£. 

In  addition,  in  the  case  of  ice  acting  on  a  vertical  or  slightly  sloping  surface 
the  crushing  strength  also  becomes  relevant.  Most  of  the  remainder 
of  this  paper  will  be  devoted  to  appraising  our  knowledge  of  these  parameters. 
For  a  more  general  review  of  the  properties  of  sea  ice  see  Pounder,  12*»13® 
Peschanskii,  19  Savel'ev, 142  Weeks,19*3  or  Zubov. 198 


2.  PHYSICAL  BASIS  FOR  INTERPRETING  SEA  ICE  STRENGTH 

Phase  Relations 


It  is  obvious  that  the  mechanical  properties  of  sea  ice  should  be  related 
to  the  amount  of  void  (brine  +  air)  in  the  ice  at  the  moment  of  testing.  There¬ 
fore,  the  first  prerequis’ for  a  detailed  analysis  of  sea  ice  data  is  an 
adequate  phase  diagram.  This  diagram  specifies  the  relative  volumes  of  ice, 
brine  and  solid  salts  that  coexist  in  sea  ice  at  equilibrium  at  a  given  tempera¬ 
ture.  The  use  of  such  a  diagram  assumes,  of  course,  that  the  ratios  of  the 
ions  in  the  system  that  is  frozen  are  the  same  as  in  standard  sea  water. 
Figure  1  shows  the  phase  diagram  that  is  currently  in  general  use  (Assur,11 
Anderson2’5).  It  gives  the  amount  of  brine,  ice,  and  solid  salt  in  the  system 
at  any  given  temperature  and  indicates  the  temperatures  of  crystallization  of 
the  solid  salts  CaC03  •  6HZ  O,  Na2  SO4  .10H2  O,  NaCl*  2H2  O,  etc.  This  dia¬ 
gram  is  based  on  the  experimental  results  of  others.132’104’105  Now  it  has 
been  known  for  some  time  that  there  were  small  but  significant  differences 
between  the  results  of  these  investigators  at  temperatures  below  -35C. 

These  differences  arc  really  not  too  surprising  considering  the  difficulties 
in  concentrating  at  these  low  temperatures  a  volume  of  brine  sufficient  to 
analyze.  These  discrepancies  were  not,  however,  considered  important 
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Figure  1.  Phase  relations  for  "standard"  sea  ice.  Circles  on  the  brine-salt 
line  indicate  temperatures  at  which  the  different  solid  salts  precipitate  (Assur). 17 


i 

i 


i 


because  the  temperature  of  natural  sea  ice  is  rarely  below  -30C.  It  was  also 
known  that  the  salts  precipitating  from  sea  ice  brine  had  never  been  positively 
identified  by  direct  (optical  or  X-ray)  methods  but  instead  had  been  deduced 
from  changes  in  the  brine  composition  and  the  general  stability  ranges  of  dif¬ 
ferent  salts  in  the  individual  pure  salt-water  systems.  In  addition,  certain 
irregularities  in  the  precipitation  of  CaCOj  *  6HZ  O  could  either  be  explained 
by  metastability  or  as  suggested  by  Tsurikov178  by  the  adsorption  of  the 
calcium  ion  on  ice.  The  possibility  of  errors  being  present  in  the  current 
phase  diagram  became  apparent  when  the  results  of  the  Russian  chemist 
Gitterman  finally  became  available  in  Savel'ev's  book.142  Savel'ev's  analysis 
of  Gitterman's  results  suggests  that  CaS04  precipitates  in  significant  quantities. 
Note  that  the  possibility  of  this  solid  salt  occurring  is  not  suggested  at  all  in 
Figure  1.  Gitterman's  work  also  indicates  that  some  of  the  early-formed  salts 
partially  rereact  with  the  brine  at  lower  temperatures  and  that  the  last  of  the 
brine  solidifies  at  -36C.  This  last  claim,  however,  appears  doubtful  inasn.ucl. 
as  Fujino52  recently  (1966)  observed  a  noticeable  change  in  the  slope  of  a  plot 
of  the  dielectric  constant  of  sea  ice  vs  temperature  at  roughly  -52C.  This 
temperature  is  close  to  -54C  where  Ringer132  suggested  that  the  last  of  the 
brine  had  solidified.  Although  there  if,  obviously  a  considerable  amount  of 
checKing  that  needs  to  be  done  on  the  sea  ice  phase  :elations,  we  are  optimistic 
in  that  we  feel  that  the  current  diagram  will  prove  to  be  reasonably  accurate  at 
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temperatures  above  the  NaCl- 2H2  O  precipitation  temperature.  This 
optimism  is  supported  by  the  recent  results  of  Richardson  and  Keller131  who 
measured  the  liquid  brine  content  directly  by  means  of  a  nuclear  magnetic 
resonance  spectrometer. 

It  would  also  be  very  important  to  establish  if  "selective  brine  drainage" 
causes  any  significant  variation  in  the  ratios  of  the  ions  in  sea  ice.  The 
current  literature141’  26’  28  on  this  subject  is  discrepant.  If  there  are  signifi¬ 
cant  changes  in  the  salt  ratios  then  the  precipitation  temperatures  of  some 
of  the  salts  will  change  and  the  applicability  of  the  standard  sea  ice  phase 
diagram  to  this  "altered"  ice  may  become  questionable. 

It  should  be  mentioned  that  there  have  been  two  direct  methods  described 
for  determining  the  amount  of  brine  in  sea  ice.  Unfortunately,  neither  the 
rather  involved  thermal  procedure  of  Savel'ev139  nor  the  nuclear  magnetic 
resonance  method  of  Richardson  and  Keller131  can  readily  be  adapted  to  field 
use. 


Volume  of  Air 


The  other  component  of  the  void  volume  is  the  volume  of  air  in  the 
sample.  It  is  relatively  simple  to  determine  the  volume  of  air  in  a  specimen 
of  sea  ice,  even  under  field  conditions. 10  »14*l49»83  Therefore,  the  dependence 
of  sea  ice  properties  on  the  air  volume  in  the  ice  was  discussed  quite  early 
in  the  Russian  literature.  Japanese  investigators159’ 148  have  also  been  able 
to  show  some  excellent  correlations  between  physical  properties  and  density 
(void  volume).  In  some  of  these  cases  it  is  quite  certain  that  the  initial  brine 
in  the  sample  had  drained  out  and  was  replaced  by  air  during  the  trip  back  to 
the  cold  room  facilities  and  the  subsequent  storage  at  relatively  high  tempera¬ 
tures  for  long  periods  of  time.  Canadian  and  American  investigators  have, 
on  the  other  hand,  preferred  to  use  the  brine  volume  computed  from  Assur's 
tables17  and  have  largely  ignored  the  volume  of  air  in  the  ice.  It  is,  of 
course,  desirable  that  both  the  volume  of  brine  and  of  air  be  determined 
independently.  Fortunately,  the  volume  of  air  is  also  quite  easy  to  determine 
for  natural  sea  ice.  The  specific  volume  of  a  sea  ice  sample  is 


where  p  is  the  density  and  A  the  relative  amount  of  the  material  specified 
by  the  subscripts  (b-brine,  s-individual  salts,  i-pure  ice,  a-air)  and  p  is 
the  density  of  sea  ice.  Because  the  amounts  of  brine  and  solid  salts  are 
proportional  to  the  ice  salinity  (weight  of  salt  per  unit  weight  of  solution) 
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b  l 
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where  the  proportionality  constants  can  be  computed  from  Assur  (Table  III). 17 
A^  can  be  represented  as 

A  =  1  -  S.  (a,  +  a  ). 
i  lbs 


(2.4) 
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Dividing  by  the  specific  volume  of  the  sample  we  find  that  the  brine  porosity 
*b'  is 


p  a.S. 
v  •  =  L-j-i 
b  Pb 

and  the  air  porosity,  v. 

c 

PA, 


(2.6) 


(2.7) 


Therefore,  the  total  porosity  v  is 


v  =  v  +  v,  1  =  1  -  p 
a  b  v 


ri  +  s.(ya.!^) 

L Pi  1  \Lf  ps  Pi  / 


(2.8) 


In  (2.8),  p  and  Si  are  measured  and  1  /p j  =  fi  (e)  and 


a,  +  a 
b  s 


=  it  (e) 


(2.9) 


can  be  calculated  as  functions  of  temperature  (0)  alone  from  Bader23  and 
Assur  (Table  HI)17  respectively.  va  can  also  be  determined  separately  by 
means  of  a  simple  experiment, 83  One  must  make  a  careful  distinction  be¬ 
tween  vjj1  in  (2.  6),  which  is  related  to  the  actual  density  p  of  the  sea  ice 
including  air  and  the  relative  brine  volume  V’D  as  introduced  by  Assur17 
which  is  related  to  the  density  of  a  theoretical  air-free  sea  ice.  The  value 
of  can  be  calculated  from  the  phase  composition  and,  with  rare  exception, 
is  the  value  used  throughout  this  paper.  In  the  future  it  would  be  desirable 
to  measure  air  porosity  separately  and  then  calculate  the  total  porosity. 

One  should  not  attempt  to  calculate  the  air  porosity  from  (2.8),  by  using 
from  tables  rather  than  v^'  for  which  tables  are  not  readily  available  a_  yet 


Structural  Considerations 


Macrostructure 


Once  the  relative  volumes  of  the  major  components  of  sea  ice  are 
established,  it  becomes  necessary  to  understand  the  structural  considera¬ 
tions  that  control  the  distribution  of  the  brine,  air  and  solid  salts  in  the 
ice.  To  do  this,  it  is  important  to  understand  how  sea  ice  forms.  In  sea 
water,  as  in  fresh  water,  the  first  crystals  to  form  are  minute  spheres  of 
pure  ice.  Growth  rapidly  changes  these  spheres  into  thin  circular  disks  in 
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Figure  Z.  Photomicrograph  of  initial  disks  and  stars  during  the  freezing  of  sea 

water  (1996,  Thule,  Greenland), 


the  general  growth  sequence  spheres  -,-disks-*hexagonal  dendritic  stars.  The 
plane  of  these  disks  or  stars  (Fig.  Z)  is  the  (9001)  or  basal  plane  which  is  the 
plane  of  maximum  reticular  atomic  density  in  an  ice  crystal.113  The  arms  of 
the  stars  form  parallel  to  the  a-axis  directions  <  1 1 Z 0>  in  the  ice  crystal 
which  is  the  shortest  lattice  vector  in  tne  (0001)  plane.  These  star-like  crys¬ 
tals  grow  rapidly  across  the  surface  of  calm  sea  water  until  they  overlap  and 
freeze  together  forming  a  continuous  thin  ice  skim.  Because  of  the  tabular 
nature  of  both  the  initial  ice  discoids  and  stars,  they  float  with  their  basal 
planes  in  the  plane  of  the  water  surface.  Ideally  this  would  produce  an  ice 
skim  with  a  completely  c-axis  vertical  orientation.  However,  since  water 
conditions  are  rarely  if  ever  completely  calm,  a  number  of  discoids  or  stars 
are  usually  caught  in  an  intermediate  position,  i.  e. ,  with  their  c-axis 
inclined  at  some  angle  to  the  vertical.  The  resulting  ice  skim  shows  a  fabric 
diagram  similar  to  that  shown  in  Figure  3:  a  slight  concentration  of  crystals 
with  their  c-axis  vertical  and  a  few  crystals  with  almost  every  other  orien¬ 
tation,  The  formation  of  surface  needles,  a  common  phenomenon  in  lake 
water,  appears  to  be  fairly  rare  in  sea  water. 

When  there  is  significant  turbulence  during  freezing,  the  effective 
thermal  supercooling  is  reduced,  more  crystals  form  per  unit  volume  of  sea 
water,  and  abrasive  action  between  crystals  is  increased.  When  the  resulting 
slush  congeals,  the  ice  cover  that  results  is  usually  several  centimeters  thick, 
fine-grained,  equigranular  in  texture,  with  a  random  c-axis  orientation.  A 


8 


THE  MECHANICAL  PROPERTIES  OF  SEA  ICE 


series  of  excellent  vertical  and  horizontal 
photomicrographs  of  such  ice  is  shown  by 
Tabata  and  Ono  (Fig.  10  and  11  -  1  to  6). 197 
This  slush  ice  should,  if  possible,  be  distin¬ 
guished  from  infiltrated  snow  ice.  The  latter 
commonly  forms  in  subarctic  regions  with 
abundant  snow  fall  when  the  weight  of  the 
overlying  snow  depresses  the  surface  of  the 
ice  sheet  below  the  water  level.  Percolation 
of  sea  water  through  cracks  changes  the  snow 
to  a  brine  saturated  slush  which  subsequently 
freezes  forming  infiltrated  snow  ice.18*  There 
currently  is  not  sufficient  information  avail¬ 
able  on  the  properties  of  slush  and  infiltrated 
snow  ice  to  know  whether  or  not  these  two  ice 
types  can  safely  be  grouped  together  when 
physical  properties  are  being  considered. 

At  some  locations  in  the  Antarctic  such  a 
considerable  portion  of  the  ice  sheet  is  com¬ 
posed  of  infiltrated  snow  ic»,  that  it  is  nec¬ 
essary  to  take  this  fact  into  account  in 
bearing  capacity  calculations. 172 

Once  a  continuous  skim  of  ice  has 
formed  across  the  sea  surface,  growth  rate 
is  determined  by  the  temperature  gradient  in  the  sheet  and  its  effective  con¬ 
ductivity.  When  a  skim  has  formed,  the  ice  crystals  lose  a  degree  of  growth 
freedom:  only  if  the  grain  boundaries  are  exactly  perpendicular  to  the 
freezing  interface  can  crystal  growth  proceed  without  one  grain  interfering 
with  the  growth  of  another.  Any  tendency  for  anisotropic  growth  will  pro¬ 
duce  geometric  selection  with  the  crystals  in  the  favored  orientation  elimi¬ 
nating  the  crystals  in  the  other  orientations  by  cutting  them  off  from  the  melt 
(Fig.  4). 114,184  Since  free-floating  ice  crystals  formed  both  from  fresh  and 
sea  water  show  pronounced  anisotropic  growth,  it  is  quite  reasonable  to 
expect  that  geometric  selection  will  occur.  A  typical  rubbir.o  of  the  bottom 
surface  of  a  Z.  5-cm  ice  skim  in  which  the  process  of  geometric  selection  is 
obviously  occurring  is  shown  in  Figure  5.  The  blank  plane-sided  polygonal 
areas  are  where  c-axis  vertical  crystals  occur.  These  crystals  grow 
appreciably  slower  than  the  surrounding  c-axis  horizontal  crystals.  They 
tl  ‘'refore  form  depressed  (~2  to  5  mm)  areas  on  the  bottom  surface  of  the 
i  skim  and  do  not  show  in  the  rubbings.  Once  the  c-axis  vertical  crystals 
e  eliminated,  rubbings  from  the  lower  levels  in  the  ice  sheet  are  similar 
to  Figure  6  which  shows  that  all  the  crystals  present  have  their  c-axes  roughly 
horizontal.  Note  that  this  process  does  not  require  the  nucleation  of  new 
crystals  with  a  c-axis  horizontal  orientation.  The  favored  orientation  forms 
by  the  survival  of  those  grains  with  this  orientation  that  are  present  in  the 
initial  ice  skim.  Bennington,  25  on  the  other  hand,  has  noted  rapid  orientation 
changes  in  the  upper  layers  of  sea  ice  sheets  that  are  difficult  to  explain  by  a 
process  of  geometric  selection.  Regardless  of  the  details  of  the  mechanisms 
involved,  it  is  clear  that  at  a  short  distance  below  the  upper  ice  surface,  sea 
ice  is  mainly  composed  of  ice  crystals  with  their  c-axes  essentially  hoiizontal. 
The  zone  in  which  this  orientation  change  occurs  has  been  called  the  transition 
zone  , 16 

Below  the  transition  zone,  s  ice  has  all  the  characteristics  associated 
with  the  so-called  columnar  zone  in  metal  ingots;181  i.e.,  a  strong  crystal 
elongation  parallel  to  the  direction  of  heat  flow',  a  pronounced  crystal  orienta¬ 
tion  (Fig.  7),  and  an  increase  in  crystal  size  over  crystals  closer  to  the  cold 


Figure  3.  Schmidt  net  plot 
(upper  hemisphere)  of  c-axis 
orientations  in  an  initial  ice 
skim  (sea  water,  31.6%c, 
Thule,  Greenland,  1956). 
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Figure  4.  Schematic  diagram  showing 
the  process  of  geometric  selection,  the 
arrows  indicate  the  direction  of  the 
c-axis  (Pe  rcy  and  Pounder).116 


Figure  5.  Rubbing  of  the  bottom  of  a 
2.5  cm  thick  ice  skim  that  is  undergoing 
geometric  selection  (Hopedale,  Labra¬ 
dor,  1955). 


Figure  fc.  Rubbing  of  the  bottom  of  a  1 -1.  8  cm  thick  sheet 
of  sea  ice;  all  c-axes  are  essentially  horizontal  (Hopedale, 
Labrador,  1955). 
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Figure  7.  Contoured  Schmidt  net  plot  of  c -axis 
orientations  at  a  depth  of  10  cm  below  the  upper 
sea  ice  surface  (Thule,  Greenland,  1957). 

Contours  represent  1,  2,  4,  6,  10  and  13  per¬ 
cent  per  1  percent  of  the  area;  200  grains. 

source.  Because  the  transition  layer  is  thin  and  occurs  at  or  near  the  ice- 
air  interface,  it  commonly  sublimes  away  or  recrysi.allizes  into  the  snow 
cover  as  the  ice  sheet  grows.  Therefore,  the  complete  thickness  of  a  one 
winter  sheet  of  sea  ice  can,  to  a  good  approximation,  be  considered  as  being 
in  the  columnar  zone. 

The  most  noticeable  change  in  the  upper  part  of  the  columnar  zone  is 
the  marked  increase  in  grain  size  with  depth  in  the  ice  sheet.  Figure  8  shows 
a  plot  of  mean  crystal  diameter  d  vs  depth  from  a  125  cm  thick  ice  sheet  at 
Thule,  Greenland.  The  measurements  were  stopped  at  60  crn  because  of  the 
large  number  of  crystals  that  were  not  completely  contained  on  a  given  thin 
section;  pronounced  linear  grain  coarsening  is  evident.  It  is  reasonable  to 
assume  that  the  crystals  stiil  present  at  the  bottom  of  a  thick  ice  sheet  not 
only  have  their  c-axes  horizontal  but  also  have  one  of  their  a-axes  close  to 
vertical. 

There  are  very  few  observations  regarding  the  grain  size  at  the  bottom 
of  thick  sea  ice  sheets.  This  lack  of  data  is,  of  course,  explained  by  the 
difficulty  in  obtaining  large  oriented  samples.  The  need  to  obtain  additional 
information  on  this  subject  is  quite  critical  in  view  of  Peyton's122’123  recent 
distinction  of  "bottom  ice"  as  a  separate  ice  type  in  the  columnar  zone. 

This  ice  is  characterized  by  a  c-axis  horizontal  orientation  and  an  extremely 
large  grain  size.  Peyton123  reports  examining  a  3  x  3  meter  block  of  1.6  m 
thick  sea  ice  and  finding  that  the  bottom  meter  exhibited  a  constant  c-axis 
orientation  over  the  entire  9  m2  cross  section.  It  also  should  be  noted  that 
in  the  old  sea  ice  incorporated  into  the  Ice  Island  Arlis  II,  Smith15'*  observed 
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Figure  8.  Mean  grain  diameter  d  vs  the 
level  below  the  upper  ice  surface  {Thule, 
Greenland,  1957). 

anisotropy  in  ice.61'*8  Highly  sutured 
intergrowths  are  common145  and  appea 
ice  grows  with  a  cellular  to  dendritic  s 


areas  as  large  as  10  m  on  a  side 
with  almost  perfect  c-axis  align¬ 
ment.  Whether  or  not  these  large 
areas  are  actually  single  crystals 
or  are  composed  of  a  number  of 
crystals  with  roughly  similar 
orientation  is  not  known. 

The  general  relation  between  the 
maximum  crystal  length  L,and  the 
maximum  crystal  width  W,  as  meas¬ 
ured  in  horizontal  thin  sections  has 
been  studied  by  Weeks  and  Hamilton188 
in  sea  ice  collected  at  Point  Barrow. 
Their  results  indicate  that  the  ratio 
L/W  is  approximately  constant  with 
a  value  of  2.  This  indicates  that  the 
rates  of  sidewise  and  edgewise 
growth  of  the  individual  crystals  do 
not  differ  nearly  as  much  as  might 
be  expected  from  studies  of  growth 
rain  boundaries  and  complex  crystal 
to  be  associated  with  the  fact  that  sea 
lid-liquid  interface. 


There  are  several  petrographic  features  of  the  columnar  zone  which  are 
best  observed  in  vertical  thin  sections.  In  contrast  to  the  parameters  dis¬ 
cussed  above,  which  are  either  constant  or  change  gradually  in  the  vertical 
(z)  direction,  these  features  commonly  appear  as  sharp  bands  which  extend 
for  large  lateral  distances.  Small  scale  horizontal  banding  in  sea  ice  has 
been  noted  in  widely  separated  regions. 158,167,85  ,25  For  example,  Shumskii 
observed  a  total  of  58  secondary  layers  in  a  291 -cm  section  of  old  pack  ice 
giving  an  average  layer  thickness  of  5  cm.  In  many  cases  these  layers  are 
produced  by  changes  in  the  amount  of  impurities  (brine  +  air)  trapped  in  the 
ice.  In  most  cases  this  layering  is  not  accompanied  by  the  nucleation  of  new 
crystals.  Similar  air  bubble  layering  in  lake  ice  is  well  known. 


The  other  type  of  horizontal  layer  is  that  produced  by  the  summer  melt 
cycle  in  multiyear  pack.  Inasmuch  as  there  are  no  observations  relating  this 
type  of  layering  to  strength  problems,  it  will  not  be  treated  here.  For  de¬ 
tailed  discussions  of  this  layering  see  Cherepanov19  and  Schwarzacner. 145 

It  is  important  that  any  determinations  of  the  physical  properties  of  sea 
ice  accurately  specify  the  type  of  ice  tested.  The  indiscriminant  mixing  of 
varied  ice  types  will  undoubtedly  make  a  significant  contribution  to  increasing 
an  already  large  natural  variability. 


Microstructure 

As  is  obvious  from  examining  a  rubbing  of  the  bottom  surface  of  a  sea  ice 
sheet  (Fig.  6),  sea  ice  crystals  possess  a  pronounced  microstructure.  The 
details  of  the  formation  of  this  substructure  are  poorly  understood  but  the  gen¬ 
eral  nature  of  the  phenomena  is  clear.  When  sea  water  freezes  the  rejection 
of  solute  by  the  solid  causes  a  solute  concentration  maximum  in  the  liquid  next 
to  the  advancing  solid-liquid  interface.  This  solute  profile,  via  the  phase 
relations,  specifies  an  equilibrium  freezing  temperature  profile  in  the  liquid. 
If  certain  temperature  profiles  occur  in  the  liquid,  a  zone  of  constitutionally 
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supercooled  liquid  is  produced  ahead  of  the  interface.  As  a  result  of  the 
formation  of  this  supercooled  zone,  a  planar  interface  ceases  to  be  stable 
and  a  cellular  or  dendritic  interface  forms.  For  a  detailed  description  of 
this  process  refer  to  references  134,  170,  171,  59,  100.  When  attempts 
are  made  to  apply  the  formal  diffusion  limited  constitutional  supercooling 
theory  to  the  natural  freezing  of  sea  and  lake  water,  severe  discrepancies 
are  noted. 1,1  These  are  undoubtedly  the  result  o±  the  fact  the  main  transfer 
mechanism  controlling  the  distribution  of  solute  in  the  liquid  ahead  of  the 
advancing  interface  is,  in  these  cases,  convection  instead  of  diffusion. 
Regardless  of  this  difficulty,  it  clearly  appears  that  the  existence  of  a  non- 
planar  interface  during  the  solidification  of  sea  water  is  the  result  of  some 
form  of  constitutional  supercooling. 

Once  a  non-planar  interface  forms,  its  exact  geometry  is  specified  by 
the  details  of  the  growth  conditions.  James72  has  discussed  the  morpholog¬ 
ical  changes  during  the  transition  to  a  fully  developed  cellular  or  dendritic 
interface.  Although  there  are  no  detailed  observations  of  sea  ice  interfaces 
during  natural  growth,  an  examination  of  the  resulting  substructure  suggests 
that  typical  morphologies  are  similar  to  the  elongated  cell  shown  in  Figure  9 
(Tiller)171  or  the  array  of  parabolic  platelets  diagrammed  by  Weeks.191  It  is 
easy  to  visualize  how  brine  becomes  trapped  in  the  intercellular  grooves  pro¬ 
ducing  the  characteristic  sea  ice  substructure  shown  in  Figures  10  and  11. 
Growth  accidents  also  may  cause  slight  orientation  differences  between  neigh¬ 
boring  plates  producing  small  angle  grain  boundaries.  Since  each  sub-grain 
represents  an  ice  platelet  that  grew  into  the  underlying  sea  water,  the 
pronounced  (0001)  elongation  of  the  plates  is  produced  by  a  slower  crystal 
growth  perpendicular  to  the  basal  plane  than  parallel  to  it.  Certain  geometric 
aspects  of  these  plates  have  been  studied  by  Weeks  and  Hamilton.188  As 
shown  in  Figure  10,  single  crystals  of  sea  ice  are  easily  delineated  under 
crossed  polaroids  because  they  behave  as  common  extinction  units.  Each 
single  crystal  is  composed  of  a  "packet”  of  plates  and  each  plate  is  separated 
by  an  array  of  brine  pockets. 

The  distance  between  adjacent  brine  layers  in  the  same  crystal  measured 
parallel  to  the  c-axis  has  been  termed  the  plate  spacing  aQ.  This  distance 
changes  systematically  with  position  in  an  ice  sheet.  That  such  changes  in 
aQ  should  be  expected  could  be  anticipated  from  the  metallurgical  literature46 
and  has  been  demonstrated  for  sea  ice,164'188  and  for  NaCl  ice.189  Figure  12 
shows  the  vertical  variations  in  the  mean  plate  spacing  5Q  for  a  sheet  of 
NaCl  ice.  The  main  factor  controlling  this  variation  is  the  general  decrease 
in  growth  velocity  as  the  ice  sheet  thickens.  Assur  and  Weeks20’  21  have 
attempted  to  fit  the  existing  NaCl  ice  data  using  a  relation  of  the  general 
form 

a0  vk  =  C  (2.  10) 

where  v  is  the  growth  velocity,  k  is  a  constant  with  a  value  of  1/2  and  C  is 
another  constant.  The  resulting  fit  is  adequate  (Fig.  13)  giving  C  =  3.  64  x 
10'4  when  a^  and  v  are  in  cm  and  cm/sec  respectively.  Unpublished  results 
of  Lofgren  and  Weeks  show  that  the  power  k  appears  to  be  a  function  of  v, 
changing  from  k~l  at  high  v  to  k~0  at  low  v.  The  use  of  k  =  1/2  does,  how¬ 
ever,  appear  to  be  a  convenient  first  approximation  over  a  wide  range  of 
growth  velocities.  More  detailed  experimental  observations  are  needed  on 
this  subject. 
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Figure  11.  Photomicrograph  of  a  thin  section  of  seaice  illustrating  brine  pocket 
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Figure  12.  Increase 
in  a0 ,  the  average 
plate  spacing,  as  a 
function  of  the  distance 
below  the  surface  of 
the  ice  sheet,  z  (Weeks 
and  Assur). 188 


Figure  13.  Logarithmic  plot  of  average  plate  spacing  a0 
vs  growth  velocity  v(Assur  and  Weeks).21  Here  (a)  is  the 
least-squares  fit,  log  a^  =  -2.9105  -  0.3729  log  v  and 
(b)  is  a  least-squares  fit  assuming  a  slope  k  of  -1/2, 
which  gives  C  =  3.  64  x  10'4  (eq  2.  10). 
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It  should  be  noted  that  currently  the  only  aspect  of  the  sea  ice  substruc¬ 
ture  that  has  definitely  been  related  to  the  growth  velocity  when  a  given 
increment  of  ice  forms  is  the  plate  spacing.  We  feci  that  it  is  very  reason¬ 
able  to  suppose  that  other  structural  features  such  as  the  systematic  spacing 
between  brine  pockets  as  measured  in  a  horizontal  plane  (see  Fig.  15)  are 
also  some  function  of  v.  As  will  be  shown  later,  it  is  very  important  that 
such  relations,  if  they  exist,  be  documented  and  incorporated  into  a  formal 
sea  ice  strength  theory. 


3.  THEORETICAL  CONSIDERATIONS 
Strength  Models 

Even  the  most  casual  examination  of  a  sea  ice  specimen  that  has  fractured 
reveals  that  the  fracture  surface  is  controlled  by  the  substructure.  Figure 
1  44  is  a  rubbing  of  a  horizontal  section  of  a  sea  ice  beam  that  failed  in 
tension.  Similar  rubbings  have  been  published  by  Tabata. 162  The  parallel 
lines  are,  of  course,  the  edges  of  individual  ice  plates.  Note  how  the  break 
only  transects  the  platelets  in  a  few  unfavorably  located  grains.  In  addition 
the  large  variations  in  measured  strength  values  for  sea  ice  made  the  develop¬ 
ment  of  formal  models  highly  desirable  in  that  they  would  suggest  compact 
forms  for  correlating  and  extrapolating  strength  results. 

The  treatment  of  the  mechanical  properties  of  sea  ice  on  the  basis  of 
structural  models  was  pioneered  by  Tsurikov, 174  In  a  subsequent  paper175 
he  modified  some  of  his  initial  calculations.  Although  he  initially  considered 
porosity  only  in  terms  of  air  content,  he  later176  introduced  brine  volume  on 
the  basis  of  Malmgren's  19^7  calculations.’7  In  his  structural  models  he 
assumed  either  uniformly  distributed  spherical  or  uninterrupted  cylindrical 
inclusions  as  voids.  Tsurikov's  imaginative  approach  also  introduced  the 
effect  of  the  orientation  of  stresses.  His  concepts  were  way  ahead  of  their 
time.  There  are,  unfortunately,  many  errors  in  principle  or  oversights 
in  Tsurikov's  work  but  his  essential  ideas  can  be  condensed  and  presented 
in  the  following  way. 


Figure  14.  Rubbing  jf  a  broken  segment  of 
sea  ice  (Anderson  and  Weeks).4 
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If  there  are  n3  spherical  voids  with  radius  r  uniformly  distributed  in 
1  cm1  of  sea  ice,  then  the  volume  porosity  is 

-  4  3  ,1 

V  -  y  IT  n3  T3  .  (3. 

It  is  essential  to  assume  that  failure  occurs  in  a  preferred  plane  which  has 
the  maximum  olane  porosity  Vp.  The  number  of  holes  in  that  plane  is 
therefore  n/'  ’  and 

v  =  ir  n3  r2  .  (3. 

P 

The  relative  strength  of  sea  ice  then  becomes 


—  =  1  -  v  n 


2/3  2 


where  <rf  is  the  failure  strength  and  <rfl  is  the  basic  strength  of  sea  ice  (see 
p.  18).  Because  r  can  be  expressed  in  terms  of  v  from  eq  3.  1 


— —  =  1  - 


where  vQ  is  the  brine  volume  for  which  the  ice  has  zero  strength. 

In  Tsurikov's  subsequent  paper,17®  he  deviated  from  this  approach  by 
assuming  that  the  solidity  of  a  line  cut  through  sea  ice  is 


where  v ^  is  the  line  porosity  (our  notation  and  definition).  This  is  in  error. 
With  a  random  distribution  of  voids,  the  line  porosity  equals  the  plane  porosity 
equals  the  volume  porosity.  In  a  systematic  array  the  corresponding  porosity 
values  depend  upon  the  arrangement.  Due  to  the  error  in  his  basic  assump¬ 
tions,  his  resulting  equations  should  not  be  used.  For  bending  Tsurikov 
reduced  the  width  and  square  of  the  height  which  is  in  error. 

For  cylindrical  voids  the  number  of  holes  (r\2  )  in  a  cut  perpendicular  to 
the  cylinders  expressed  per  cm2  becomes  important.  The  plane  porosity 
equal  to  the  volume  porosity  is 


-  v  - 


(3.7) 


THE  MECHANICAL  PROPERTIES  OF  SEA  ICE 


17 


Again  assuming  failure  in  a  preferred  plane  parallel  to  the  cylinders,  the 
maximum  line  porosity  v f  becomes  important.  The  number  of  holes  in  that 
line  is  ‘•Jn^  ar.d  the  line  porosity  is 

vj  =  2  \Jnl  r.  (3.8) 


The  relative  strength  of  sea  ice  in  a  plane  parallel  to  the  cylinder  becomes 


—  =1-2  \l  n2  r 

(ro 


and  substituting  for  r  from  eq  3.  7 


<r 

o 


with 


..  IT 

vo  ~  4  • 


(3.9) 


(3.  10) 


(3. 11) 


The  corresponding  equation  in  Tsurikov174  is  wrong  by  oversight  but  was 
corrected  later.175  For  failure  in  a  plane  perpendicular  to  the  cylinders, 
the  strength  is  simply 


(3.12) 


which  is  always  higher  than  eq  3.  10. 

Other  Russian  investigators  appear  to  have  paid  little  attention  to 
Tsurikov's  work  although  it  has  been  referenced  in  a  recent  Russian  review70 
of  American  work  in  sea  ice  physics.  Tsurikov's  papers  clearly  indicate  the 
possibility  of  developing  a  structurally  based  theory  of  sea  ice  strength.  The 
next  step  one  must  take  is  a  proper  consideration  of  the  platelet-like  micro¬ 
structure  of  sea  ice  crystals  where  brine  pockets  are  arranged  in  systematic 
parallel  arrays. 

This  step  was  provided  by  Assur, 17  Anderson  and  Weeks,4  and 
Anderson2’5  who  fortunately  had  available  a  considerable  amount  of  informa¬ 
tion  on  die  structure  of  sea  ice.  They  noted  that  in  a  relatively  cold  sample 
of  sea  ice  (Fig.  15)  with  a  small  brine  volume,  the  brine  pockets  are  small 
and  widely  spaced.  When  a  potential  failure  plane  between  two  ice  plates  is 
considered,  it  is  found  that  an  appreciable  percentage  of  this  plane  is  ice. 
Therefore,  this  sea  ice  would  be  expected  to  have  a  significant  tensile 
strength.  The  specimen  shown  in  Figure  16,  on  the  other  hand,  contained 
a  very  large  volume  of  brine  while  it  was  located  in  the  ice  sheet.  This 
brine  was  later  lost  during  storage  and  preparation.  A  plane  can  be  passed 
through  this  sample  without  encountering  more  than  a  few  ice-ice  bonds. 


Figure  15.  Photomicrograph  of  sea  ice  at  low  temperatures  (Thule,  Greenland, 

1957). 


Therefore,  this  specimen  would  be  expected  to  have  a  tensile  strength  close 
to  zero  even  though  it  contains  an  appreciable  amount  of  ice.  This  is 
exactly  what  is  observed  in  strength  tests  (see  Fig.  ZZ):  the  failure  strength 
orf  decreases  from  a  maximum  <r0  at  zero  brine  volume  (vj,  =  0)  to  a  strength 
of  zero  at  a  brine  volume  vQ  .  Therefore,  it  should  be  possible  to  express 
sea  ice  strength  in  the  general  form 


where  crQ  can  be  considered  the  basic  strength  of  sea  ice  (i.  e.  ,  the  strength 
of  an  imaginary  material  that  contains  no  brir.e,  but  still  possesses  the  sea 
ice  substructure  and  fails  as  a  result  of  the  same  mechanism  that  causes 
failure  in  natural  sea  ice)  and  ij;  is  the  "plane  porosity"  or  relative  reduction 
in  area  of  the  failure  plane  as  the  result  of  the  presence  of  brine  and  air 
inclusions.  The  critical  value  of  4*  in  the  failure  plane  is 


where  v  is  the  void  volume  or  porosity  and  va  and  are  the  volume  of  air 
and  brine  respectively  in  the  ice.  We  will  now  for  simplicity  consider  a  sea 
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Figure  16.  Photomicrograph  of  sea  ice  at  a  temperature  only  slightly  below  the 
freezing  temperature  (-3C,  Thule,  Greenland). 

ice  specimen  in  which  »  ^a.  To  express  4<  in  terms  of  vjj,  a  simplified 
model  of  the  geometry  of  the  brine  pockets  must  be  given.  This  can  be 
done17  as  shown  in  Figure  17.  Here  the  relative  brine  volume  is 


b  a0  b0  g0 

where  the  symbols  are  defined  in  Figure  17.  Now  defining 


(3. 17) 


we  find  that  the  reduction  in  cross  sectional  area  as  the  result  of  the 
presence  of  the  brine  pockets  is 

=  2rbg  =  2rbV 
8e  Po  ao 


(3.  18) 


Figure  17.  Idealized  diagram  of  the  shape  of  the 
brine  inclusions  in  sea  and  NaCl  ice  (Assur).11 


The  effective  strength  due  to  this  reduction  in  cross  sectional  area  becomes 


°J_  =  i  2rb  V 
o'o  Po  ao 


(3.  19) 


The  question  now  is  how  the  geometric  parameters  r^  and  y  vary  with 
Here  a  number  of  different  assumptions  can  be  made.  The  simplest  case 
is  to  assume  that  geometric  similarity  is  preserved  only  along  the  B-axis 
(see  Fig.  17)  and  that  the  width  a3  well  as  the  length  of  the  brine  pockets 
remains  constant.  In  this  case  rj,  must  change  proportionally  to  V5  and  an 
equation  of  the  form 


o-o 


(3.20) 


results.  If,  on  the  other  hand,  the  average  length  and  spacing  of  the  brine 
cylinders  remain  constant,  then  changes  in  t/j,  will  be  reflected  only  in  the 
BC  cross  section.  If  geomotr  similarity  is  preserved  in  this  cross 
section,  then  our  equation  is  of  the  form 


h 

<T0 


1/2 


(3.21) 


If  all  brine  pockets  remain  of  a  similar  shape  during  changes  in  v^,  all 
linear  dimensions  will  change  proportionally  to  (geometric  similarity 

in  space).  This  results  in  an  equation  of  the  form 


-  c 


2/3 
vb  • 


1 


(3.2*) 
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These  three  sets  of  models  can  be  represented  as  straight  lines  in  o-f,  vk 
coordinates  where  k  is  1,  1/2  and  2/3  respectively.  The  o-f  axis  intercept 
is  ng  and 


(3.23) 


where  v0  is  the  volume  of  brine  necessary  to  cause  the  ice  to  have  zero 
strength.  Defining 


(3.24) 


(3.25) 


vt>  for  the  three  classes  of  models  may  be  written 
.  _  -  _  1  _  F 

k  =  1,  V„  -  -  -  - ~ 

0  P<10  3  r  a0 


(3.26) 


(3.27) 


(3.28) 


A  large  number  of  specific  models  can  be  developed  under  these  three 
categories17''4.  The  two  that  have  been  utilized  in  discussing  specific 
strength  results  are  t.he  constant  width  and  the  elliptical  cylinder  models. 

In  the  constant  width  model 


F  =  4  r  ru 
a  b 


(3.29) 


J.  =  d^ 


(3.30) 


where  d0  is  the  minimum  width  of  a  parallel  brine  layer  before  it  splits  as 
the  result  of  interfacial  tension  to  produce  individual  brine  pockets  It  is 
usually  assumed1®9’ 4  that  d^  ~constant  because  the  ice-brine  interlacial 
tension  would  not  be  expected  to  be  particularly  temperature-sensitive  in 
the  temperature  interval  where  sea  ice  observations  are  made.  A  value 
of  d0  of  7  x  10*3  cm  has  been  measured  in  natural  sea  ice4  by  direct 
petrographic  methods  and  a  value  of  1 1  x  10"3  cm  in  NaCl  ice  by  an  indirect 
calculation. 189  When  a  brine  layer  becomes  thinner  than  this  value  it  i'.ecks, 
causing  ice-ice  bonds  to  form.  When  the  lirst  of  these  bonds  forms  between 
the  growing  plates  at  the  bottom  of  an  ice  sheet,  it  establishes  the  location 
of  the  so-called  bridging  layer.  The  ice  platelets  below  this  location  are 


mam 
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not  laterally  connected  and  therefore  have  zero  tensile  strength.  This 
strengthless  layer  has  been  termed  the  skeleton  layer  since  it  can  bt 
considered  to  be  the  result  of  "skeletal"  growth. 


In  the  elliptical  cylinder  model  let 


(3.31) 


ir  rr 
b  a 


(3.  32) 


and  eq  3.21  becomes 


—  =1-2  J-LX  -JV7 

(To  '  tPo  b 


(3.33) 


If  we  wish  to  consider  circular  cylinders  (e  =  1)  and  neglecting  any 
interruption  of  brine  pockets  in  the  vertical  direction  (y  =  1),  eq  3.  33 
becomes 


0o  1 


•'/"  Po 


(3. 34) 


As  will  be  discussed  in  detail  later,  recent  experimental  results50 
have  verified  the  tentative  suggestion  of  Assur17  that  the  strength  of  warm 
sea  ice  does  not  approach  zero  as  a  linear  function  of  v^  as  suggested 

rhese  previous  models.  Instead,  if  the  test  results  are  plotted  on  o-f, 
coordinates,  the  resulting  curves  appear  to  be  convex  toward  the 
origin  of  the  coordinate  system  in  the  high  brine  part  of  the  diagram  (Fig. 
26).  It  is  possible  to  modify  eq  3.  33  in  such  a  way  that  it  more  properly 
describes  this  observed  strength  variation 


L  .  .  -  ills 


(3. 35) 


where 


(3.  36) 


with  a  value  of  *7 v0  ~  0.  5  if  *7 v0  is  computed  from  the  linear  intercept. 
We  may  choose  a  power  n  >  1  for 


0o  1 


(3.  37) 
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without  losing  the  advantages  of  the  theoretical  models.  The  simplest 
possible  assumption  is  n  =  2  which  gives  an  expression  that  can  be 
linearized  in  the  form 


(3.38) 


where  by  plotting  \f<rf  vs  n/v^  ,  we  obtain  \/^ as  the  intercept  and  \/aro  7vq 
as  the  slope.  From  eq  3.  37  with  n  =  2  we  obtain 


which  for  =  v0  equals  zero.  Therefore,  the  proper  choice  of  v0  may  not 
be  important  since  the  change  in  <r£  in  its  vicinity  is  small.  Equation  3.  39 
expresses  a  gradual  vanishing  of  <rf  rather  than  an  abrupt  decrease.  Equation 
3.  37  with  n  -  2  can  also  be  written  in  the  form 


<r0 


(3.40) 


which  is  a  second  order  polynomial  in  v/v^  .  Here  we  have  a  combination 
of  the  two  main  models  discussed  above;  one  is  linear  in  s/v^ ,  the  other  in 
vjj.  F rom  eq  3.  40 


d(<If/  <r0) 
dv/vo 


2 

-/vo 


v  o  b 


(3.41) 


where  for  =  0  exactly  and  for  small  as  an  approximation  gives 


d(<rf/a-0)  2 

ds/vo 

with  the  simplest  possible  assumption  v0  -  1,  we  have 


^(«rf/  (ro  ) 

- 1 -  =  -  2 

dN/v^ 

which  is  in  very  good  agreement  with  observation. 17,187,50 


(3.42) 


(3.43) 


Stress  Concentration  and  Brittle  Fracture 

By  reading  Assur,  17  Anderson  and  Weeks,  4  Anderson,  2,5  and  the 
discussion  in  Natioral  Academy  of  Sciences  -  National  Research  Council 
598,  it  becomes  apparent  that  there  has  been  considerable  diversity  of 
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opinion  regarding  certain  aspects  of  the  problem  of  stress  concentration  in 
sea  ice.  The  main  portion  of  this  discussion  was  involved  with  Anderson 
and  Weeks4  presenting  eq  3.13  in  the  form 

°f  =  TT  (l  '  W  (3- 44* 

where  is  the  strength  of  pure  ice  and  k  is  a  stress  concentration  factor 
caused  by  the  presence  of  brine  pockets  in  sea  ice.  For  small  isolated 
circular  brine  pockets,  k  was  assumed  to  be  ~  3.  This  implied  that  the 
basic  strength  of  sea  ice  g0  was  1/3  the  strength  of  lake  ice  and  that  this 
difference  was  produced  by  the  presence  of  brine  pockets  in  sea  ice.  That 
this  was  incorrect  was  shown  experimentally  when  Weeks187  and  Graystone 
and  Langleben57  demonstrated  that  ~  <j-q  for  NaCl  ice  and  sea  ice  respec- 
cively.  This  does  not  mean  that  stress  concentrations  are  not  produced  by 
the  sea  ice  microstructure  but  only  as  suggested  by  Assur17  ,  that  similar 
stress  risers  also  occur  in  fresh  water  ice.  The  factor  of  ~  3  observed 
by  Anderson  and  Weeks4  was  produced  by  differences  in  stress  concentra¬ 
tion  associated  with  differences  in  testing  procedures  and  not  with  differences 
in  k  associated  with  the  presence  or  absence  of  the  microstructure.  In 
actuality,  the  equating  of  with  <j-j  is  a  suspect  procedure.  53  Inasmuch  as 
the  microstructure  of  the  fresh  water  ice  is  unspecified,  crj  is  undetermined. 
In  addition  fresh  water  ice  which  by  definition  contains  no  brine  pockets  quite 
possibly  fails  by  a  different  mechanism  than  does  sea  ice.  54 

Goetze53  has  pointed  out  that  it  also  is  possible  to  approach  the  sea  ice 
strength  problem  as  a  problem  in  elastic  instability  in  a  brittle  material. 

The  well  known  Griffith  crack  problem  is  of  this  type.  In  the  Griffith  crack 
analysis,  it  has  been  shown  that  two  conditions  must  be  satisfied  for  sponta¬ 
neous  crack  growth: 

1.  A  stress  concentration  must  exist  that  is  sufficient  to  overcome  the 
intrinsic  strength  of  the  material  and 

2.  Elastic  readjustment  as  the  crack  develops  must. provide  a  de¬ 
creasing  energy  path  for  the  specimen  as  a  whole. 

Griffith's  criterion  for  a  failure  stress  is 


where  erf  is  the  failure  stress,  E  is  Young's  modulus,  y'  is  the  total  energy 
required  per  unit  increase  in  the  area  of  the  crack  (comprised  of  the  true 
surface  energy  and  the  plastic  work  done  by  the  stress  concentration  at  the 
tip  of  the  moving  crack),  and  c  is  the  crack  length.  Goetze  has  extended 
this  treatment  by  considering  the  failure  criterion  for  a  polycrystalline, 
elastic  material  which  contains  a  "microstructure"  capable  of  generating 
a  stress  concentration.  This  is  done  by  considering  a  grain  that  has  for 
some  "eason,  perhaps  orientation  differences,  an  apparent  E  lower  than 
its  surroundings.  The  result  is 
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where  G  is  the  shear  modulus,  t)  =  ( 1  -  p.)  /  ( 1  +  p.)  where  p  is  Poisson's  ratio, 
K  is  the  proportion  of  the  stress  <r  that  is  exerted  on  the  weak  crystal 
( 0  <K <_  I),  and  6  is  a  complicated  function  of  c0  and  c  where  2c0  is  the 
original  grain  length  and  (c  -  c3)  is  the  length  of  the  crack  (prior  to  cracking 
c  =  c0  and  6  =  it).  If  (1  -  K)  >  0,  failure  can  clearly  occur  for  sufficiently 
large  <r  without  the  presence  of  a  crack  prior  to  failure.  This  is  not  true  in 
a  homogeneous  material.  Goetze  then  shows  by  a  dimensionless  argument 
that  elastic  instability  can  occur  during  plane  stress  when 


<r  > 


V 


r^v7 


(3.47) 


where  P  has  a  value  close  to  1  and  when  a  class  of  materials  is  considered 
in  which  the  microstructure  varies  only  in  size,  P  remains  constant.  The 
length  d  is  taken  to  be  a  characteristic  of  the  size  of  the  microstructure. 
Therefore,  when  dealing  with  failure  stresses 
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(3.48) 


Considering  ring  tensile  results  on  sea  ice,  we  find  that  oy  ~  1 5  kg/cm2 
and  that  E  ~5  -  10  x  1010  dynes/cm2  ,  y'  ~76  ergs/cm2  and,  assuming  that 
P  ~1,  eq  3.48  implies  that  the  fracture  originates  in  a  microstructure  with 
a  size  of  0.  1  to  1.0  mm.  This  is  clearly  of  the  same  scale  as  the  brine 
pockets  in  sea  ice. 

To  treat  sea  ice  as  ar.  instability  problem  it  is  necessary  to  assume 
that  only  the  brine  pockets  within  a  given  plane  interact  elastically.  This, 
of  course,  is  quite  reasonable  considering  the  structure  of  sea  ice.  We 
then  must  obtain  the  elastic  solution  for  a  periodic  array  of  cracks.  This 
apparently  is  not  known.  Therefore,  Goetze  makes  a  number  of  approxima¬ 
tions  to  get  some  idea  of  the  general  solution.  If  bp  and  r^are  defined  as 
shown  in  Figure  17,  his  approximate  fracture  relation  is 

r2"  -  * 

2rb 

He  then  assumes  that  2  r^  can  be  identified  with  n/vl  where  is  the  brine 
volume  and  shows  that  a  oy  vs  1  /'Jv^  plot  of  Assur's17  ring  tensile  data  is 
to  a  good  approximation  linear.  There  are  some  obvious  difficulties  with 
eq  3.49  -  note  that  when  r^  —0,  ay  —  w.  Yet  it  is  cold  sea  ice  with  its 
correspondingly  low  brine  volumes  that  most  closely  approximates  the 
brittle  behavior  postulated  by  an  elastic  instability  analysis.  Nevertheless 
this  is  an  interesting  and  promising  approach,  if  not  for  sea  ice  which  is 
always  quite  close  to  its  melting  temperature,  then  for  other  porous 
materials  which  show  brittle  fracture  characteristics.  Work  along  these 
lines  should  be  continued. 


ay  >,  const 


(3.49) 


Air  Bubbles  and  Salt  Reinforcement 


It  should  be  noted  that  additional  petrographic  observations  are  needed 
on  the  distribution  of  air  pockets  in  sea  ice.  If,  as  seems  reasonable,  air 
bubbles  are  localized  along  the  same  substructure  as  the  brine  pockets,  then 
we  may  simply  treat  va  as  an  additional  contribution  to  and  deal  with  the 
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total  void  volume  v  where  v  =  va  +  vjj.  If,  on  the  other  hand,  the  geo¬ 
metrical  considerations  governing  the  distribution  of  the  air  bubbles  are 
quite  different  from  those  controlling  the  brine  distribution,  va  will  have 
to  be  considered  separately. 

There  has  been  some  discussion  in  the  recent  sea  ice  literature  re¬ 
garding  the  possible  effects  of  Na2  SO4  •  I0H2  O  and  NaCI*  2H2  O  in  increasing 
the  values  of  a- f  beyond  the  normal  increase  that  would  be  associated  with 
the  decrease  in  the  size  of  the  brine  pockets  as  the  solid  salt  precipitates. 
Possible  hypotheses  for  this  so-called  solid  salt  reinforcement  of  the  brine 
pockets  have  been  advanced  by  Assur17  and  Peyton.123  Because  the  experi¬ 
mental  verification  for  this  effect  is  still  in  doubt,  we  will  not  discuss 
these  theories  here.  It  shjuld,  however,  be  noted  that  both  hypotheses 
depend  upon  how  the  solid  salts  are  distributed  in  the  brine  pockets  during 
precipitation.  The  salts  may  precipitate  in  a  solid  mixture  with  the  ice 
in  which  case  reinforcement  is  possible.  Or  the  solid  salt  crystals  may 
remain  in  the  liquid  brine  in  which  case  no  reinforcement  is  possible.  If 
the  solid  salts  enter  into  a  solid  mixture  with  ice,  they  will  not  participate 
in  further  brine  drainage  and  a  change  in  the  ion  ratios  with  time  should 
be  observed.  If  the  solid  salts  drain  together  with  the  brine,  the  ion  ratios 
should  remain  constant  in  sea  ice.  Sporadic  evidence  to  support  either  of 
these  views  is  available.26,28, 142 


Interrelations  Between  Growth  Conditions  and  Strength 

Assur  and  Weeks20’21  have  discussed  how  changes  in  the  geometry  of 
the  sea  ice  microstructure  as  produced  by  changes  in  the  meteorological 
conditions  during  growth  can  influence  the  failure  strength.  Tc  do  this 
one  may  assume  that  eq  2.  10,  the  relation  between  the  plate  spacing  a0  and 
the  growth  velocity  (Fig.  13),  can  be  expressed  as 

a0  sJV  =  2  \(cD  Az  ~  const.  (3.50) 

Here  D  is  the  molecular  diffusion  coefficient  of  salt  in  water,  c  is  a 
coefficient  which  accounts  for  turbulent  transfer  of  solute  in  excess  of 
molecular  diffusion  and  Az  is  a  characteristic  length  associated  with  the 
distance  a  platelet  grows  before  the  salt  diffusion  for  neighboring  platelets 
meets  between  the  platelets  preventing  another  platelet  from  protruding. 

It  should  be  noted  that  the  fact  that  details  of  the  justification  for  eq  3.  50 
are  undoubtedly  incorrect  is  not  important  in  the  following  discussion  as 
long  as  some  relation  of  this  general  form  holds.  Now  the  fresh  ice  growth 
equation  is 
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where  t  is  the  time,  h  the  ice  thickness,  p  the  density  of  ice,  L  the  latent 
heat  of  fusion  of  ice,  A  0  the  difference  between  the  ambient  temperature  and 
the  freezing  temperature  of  the  water,  e  the  overall  coefficient  of  surface 
heat  transfer,  and  K  the  thermal  conductivity.  This  relation  may  be  applied 
to  sea  ice  if  slight  adjustments  are  made  in  the  values  of  L  and  K.  Now 
rearranging  and  differentiating  eq  3.51 
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which  when  combined  with  eq  3.  50  gives 
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(3.53) 


For  thin  ice  1/e  »  h/K  and  a0  is  independent  of  the  ice  thickness  h,  while 
for  thick  ice  h/K»  1/e  giving  a0  varying  proportional  to  ^/K.  The 
empirical  power  of  1/6  as  determined  by  Weeks  and  Assuf,w  lies  between 
these  theoretical  limits  of  zero  and  1/2  (Fig.  12). 

Because  a^  appears  in  our  model  strength  equations  we  may  substitute 
in  3.  20  and  3.  30  giving 
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which  for  thin  ice  gives 


’Ll 

Oo 


const 


vb 

n/A0 


and  for  thick  ice 


(3.54) 


(3.55) 
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where  z  is  the  distance  of  a  given  ice  layer  below  the  upper  ice  surface. 
Equation  3.55  provides  a  scientific  justification  for  the  intuitive  conclusions 
of  early  Russian  observers  who  felt  that  the  strength  of  young  sea  ice  was 
in  some  way  influenced  by  its  growth  conditions.  Assur  and  Weeks20,21  also 
consider  the  effect  of  a  snow  cover  and  ways  of  incorporating  the  change  in 
the  salinity  of  the  ice  sheet  with  time.  Although  it  is  beyond  the  scope  of 
this  review  to  discuss  the  sea  ice  salinity  profile  and  its  time  dependence, 
an  appraisal  of  the  current  knowledge  of  this  subject  can  be  obtained  from 
references  178,  28,  190  and  179. 


4.  EXPERIMENTAL  RESULTS 

The  quality  of  experimental  studies  on  the  mechanical  properties  of 
sea  ice  is  quite  varied.  The  fact  that  a  large  number  of  published  papers 
do  not  contain  sufficient  supplementary  information  to  permit  an  analysis 
of  strength  results  was  noted  by  Zubov  in  1945  and  unfortunately  is  still 
true.  Ideally  every  measurement  of  a  mechanical  property  should  at  least 
be  coupled  with  a  measurement  of  the  salinity  Sj  and  temperature  Qi  of  the 
tested  sample  (v^  =  f(Sj0j)).  A  measurement  of  the  air  temperature  9ais 
not  adequate.  It  is  very  perplexing  to  attempt  to  analyze  data  that  are 
accompanied  only  by  a  statement  that  0a  =  +4C.  In  addition,  the  orientation 
of  the  sample,  the  rate  of  stress  application  a.  the  volume  of  included  air 
va,  details  of  the  exact  location  of  the  sample  in  the  ice  sheet,  the  history 
of  the  ice  sheet  and  the  handling  of  the  sample  after  it  is  removed  from  the 
sheet  should  be  noted.  It  is  also  desirable  to  have  at  least  a  limited 
petrographic  description  of  the  structural  aspects  of  the  sample.  In  this 
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review  we  will,  in  general,  only  discuss  those  papers  which  have  provided 
reasonable  supplementary  information.  Exceptions  will  only  be  made  in 
cases  such  as  shear  strength  where  there  are  almost  no  data  available  re¬ 
gardless  of  quality.  We  have,  however,  in  the  references  attempted  to 
give  a  fairly  complete  list  of  the  papers  that  consider  any  aspect  of  the 
mechanical  properties  of  sea  ice. 

General  descriptions  of  recommended  test  procedures  and  equipment 
that  has  been  used  in  studying  the  mechanical  properties  of  sea  ice  in 
Russia  and  in  the  United  States  and  Canada  are  given  in  Kudriavtsev82  and 
Peschanskii,119  and  Butkovich35  respectively.  The  main  types  of  tests 
considered  are  tensile  (ring  and  "dog-bone"),  flexural  (small  beam  and  in- 
situ)  and  shear,  Kudriavtsev  also  discusses  an  indentation  test  and  a  rather 
involved  determination  of  the  shear  modulus  and  the  coefficient  of  viscosity 
by  torsion.  These  two  tests  have  apparently  not  been  used  on  sea  ice.  Al¬ 
though  there  are  some  differences,  the  test  procedures  suggested  are,  in 
general,  similar.  Particular  stress  is  placed  on  the  importance  of  de¬ 
termining  the  salinity  and  temperature  of  the  sample,  ar.d  it  is  recommended 
that  for  tests  concerned  with  the  elastic  behavior  of  sea  ice  the  rate  of  stress 
application  <r  should  be  greater  than  0.5  kg/cm2  -  sec.  This  particular  value 
of  <r  is  based35  on  the  observation  by  Jellinek73  that  the  tensile  strength  of 
snow  ice  cylinders  is  independent  of  a  at  o->  0.5  kg/cm2  -  sec  (Fig.  18). 
Evidence  for  pure  ice  also  shows  this  to  be  a  reasonable  value. 

The  importance  of  testing  sea  ice  immediately  after  it  has  been  removed 
from  the  ice  sheet  should  be  stressed.  Exceptions  to  this  procedure  should 
only  be  made  when  storage  facilities  with  ambient  temperatures  of  -30C  or 
colder  are  available.  If  the  sample  is  stored,  the  coaxial  hole  should  not 
be  drilled  until  shortly  before  testing.  If  these  procedures  are  not  followed 
brine  drainage  invariably  results  in  the  surface  layer  of  the  specimen  where 
failure  is  initiated.  If  samples  are  stored  for  appreciable  lengths  of  time  at 
temperatures  warmer  than  -23C,  extensive  brine  drainage  throughout  the 
complete  sample  results.  When  this  occurs,  the  correlation  between  the 
mechanical  properties  of  the  ice  and  its  brine  volume  disappears  and  is  re¬ 
placed  by  a  "fictitious"  correlation  between  mechanical  properties  and 
sample  density  as  observed,  for  instance,  by  Tabata.  15*’1M 

If  a  test  temperature  other  than  the  in-situ  temperature  is  desired,  it 
is  necessary  to  let  the  specimen  sit  for  approximately  3  hours  before  the 
temperature  change  at  the  center  of  the  core  is  95%  complete.126  Conversely 


,  «Tl  Of  STHESS  (kg/cm'-mci 

Figure  18.  Average  tensile  strength  as  a  function  of  rate  of 
stress  application  !r  for  snow  ice  cylinders;  height  2  cm,  diam 

2  cm  (Jellinek.)73 
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if  an  in-situ  test  temperature  is  desired  and  it  is  necessary  to  test  at  a 
different  temperature,  the  center  temperature  of  the  core  can  be  assumed 
to  be  within  5%  of  the  original  temperature  if  the  test  can  be  performed 
within  approximately  10  minutes  of  the  time  of  removal  of  the  specimen. 

No  information  is  available  on  the  rate  of  change  of  the  ice  surface  tempera 
ture  which  is  an  important  factor  in  strength  tests. 


Tensile  Strength 


Ring  tensile  tests 

The  ring  tensile  test  has  been  extensively  applied  to  studies  of  sea  ice. 
There  are  several  reasons  for  this:  the  SIPRE  coring  auger  produces  a 
suitable  sample,  the  test  is  simple,  easy  to  perform,  and  shows  pronounced 
changes  with  temperature,  salinity,  depth  and  microstructure.  The  teBtwas 
initially  designed  to  study  the  tensile  strength  of  such  brittle  materials  as 
concrete  and  rock.  The  theory  for  the  test  was  developed  in  an  article  by 
Ripperger  and  Davids113  and  further  elaborated  by  Assur. 17  The  test 
specimen  consists  of  a  sea  ice  cylinder  with  a  3-in.  (7. 6-cm)  diam  in  which 
a  1/2-in.  ( 1 . 27 -cm)  to  1  -in.  (2.  54-cm)  coaxial  hole  has  been  drilled.  The 
specimen  is  then  subjected  to  a  compressive  force  normal  to  the  axis  of  the 
cylinder.  The  sample  fails  as  a  result  of  tensile  stresses  that  develop  in 
a  vertical  plane  with  a  maximum  at  the  inner  hole  perpendicular  to  the 
direction  of  the  load.  The  equation  for  computing  the  ring  tensile  strength 
is 
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KP 
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where  P  is  the  load  at  failure,  f  is  the  average  length  of  the  sample,  and 
rQ  is  the  outer  radius  of  the  sample.  K  is  a  stress  concentration  factor 
that  is  a  complicated  function  of  the  ratio  of  the  inner  and  outer  radii 
( r  i /  r0)  °f  hollow  cylinder.  The  theory  for  the  test  predicts  K  =  1  for 
a  solid  cylinder  without  a  hole.  For  an  infinitesimally  small  hole  K  =  6 
and  for  rj/rQ  =  1/6,  which  has  commonly  been  used,  K  =  7.09.  Studies 
of  the  effect  of  variation  of  the  hole  size  on  the  resulting  strength  have 
been  made  by  Butkovich37  and  Frankenstein.1*0  Butkovich  found  no  change 
in  the  values  of  <r£  for  glacier  icc  when  the  ratio  (rj/rQ)  was  changed  from 
1/3  to  1/6.  Frankenstein  compared  sea  ice  tests  performed  on  specimens 
with  a  radial  hole  diameter  of  1/2  in.  (1.27  cm)  with  tests  performed  on 
solid  sea  ice  cylinders  without  a  predrilled  hole  but  containing  brine  pock¬ 
ets.  Assuming  that  the  solid  cylinder  should  give  the  same  o-f  value  as  the 
predrilled  specimen,  a  new  K  with  a  value  of  5.  2  was  computed.  Ice 
temperature  anu  brine  volume  do  not  appear  to  affect  this  value. 

Because  of  the  geometry  of  the  ring  tensile  test,  initial  failure  is 
restricted  to  a  very  small  volume.  It  is  difficult  to  estimate  this  volume 
exactly,  but  an  approximate  calculation53  suggests  a  region  1.5  mm  x 
0.  3  mm  x  the  length  of  the  cylinder.  In  short,  ring  tensile  tests  force 
failure  upon  a  volume  that  has  dimensions  only  slightly  larger  than  the 
dimensions  of  the  sea  ice  microstructure. 


The  theory  predicts  sharp  stress  concentrations  in  this  small  volume. 
No  doubt,  stress  relief  occurs  as  a  result  of  plastic  deformation  which 
should  be  a  function  of  temperature.  This  causes  exaggerated  values  and 
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additional  difficulties  in  interpretation.  Sala137  has  discussed  this  difficulty 
but  has  not  contributed  sufficient  experimental  evidence  to  document  his 
suggestion  that  the  stress  concentration  index  for  natural  sea  ice 


= 


(4.2) 


where  a  is  the  theoretical  stress  concentration  factor  and  (3  the  reduced 
stress  concentration  factor  as  observed  in  experiments,  is  close  to  zero. 

Butkovich14  was  the  first  investigator  to  apply  the  ring  tensile  test  to 
the  study  of  sea  ice  at  Hopedale,  Labrador.  He  found  that  the  fai'.ure 
strength  <rf  did  not  appear  to  be  a  strong  function  of  either  the  sample 
density  or  its  position  in  the  ice  sheet.  He  also  showed  that  a,  increased 
with  decreasing  ice  temperature  ^  in  the  temperature  range  -2,  5  to  -19.  1C 
for  sea  ice  specimens  with  comparable  densities  and  salinities.  In  a  later 
series  of  tests  at  Thule.  Greenland,  Butkovich36  further  documented  the 
increase  in  07  with  a  decrease  in  .  He  also  showed  that  af  decreased 
as  the  brine  volume  vf,  increased  (Fig.  19).  Samples  obtained  by  horizontal 
coring  which  were  oriented  so  that  failure  occurred  across  the  micro- 
structure  gave  consistently  higher  af  values  than  values  obtained  from  verti¬ 
cal  cores  with  no  specific  orientation  to  the  microstructure.  As  Butkovich 
pointed  out,  this  difference  in  af  is  quite  reasonable  in  terms  of  the  struc¬ 
ture  of  the  sea  ice.  He  also  noted  that  for  identical  values  of  vj,  and 
Thule  ice  gave  consistently  higher  <rf  values  than  Hopedale  ice.  A  possible 
explanation  was  advanced  for  this  difference:  that  af  is  in  some  way  influ¬ 
enced  by  the  past  thermal  history  of  the  sample.  A  similar  hypothesis  had 
also  been  advanced  by  Assur17  who  proposed  that  strength  changes  might 
lag  appreciably  behind  changes  in  Gj.  Attempts  to  resolve  these  suggestions 
by  studying  the  Thule  data  were  inconclusive. 
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Figure  19.  Ring  tensile  strength  of  sea  ice  vs  brine  volume,  vertical  cores. 
Circled  points  are  group  averages  of  10  or  more  tests  (Butkovich).36 
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Assur17  reanalyzed  some  of  Butkovich's  data  in  addition  to  new  data, 
in  terms  of  a  o-f  vs  \Zvjj  plot  which  was  suggested  by  the  development  of 
strength  models  as  discussed  earlier.  His  final  graph  is  shown  as  Figure 
20.  From  this  graph  and  accompanying  theoretical  arguments,  he  advanced 
a  number  of  hypotheses  regarding  the  strength  behavior  of  sea  ice.  The 
more  important  of  these  were 

1.  At  high  temperatures  (above  -8.  2C),  the  strength  of  sea  ice  signifi¬ 
cantly  decreases  with  increase  in  (Fig.  20,  line  Bf-B2  ). 

2.  At  temperatures  between  -8.2  and  -22. 9C,  the  presence  of  solid 
Na2  SO4 '  10H2  O  increases  the  strength  of  the  ice  by  roughly  one  third 
(Fig.  20,  line  Cj-C). 

3.  At  temperatures  below  -22. 9C,  the  precipitation  of  NaCl*  2H2  O 
causes  sea  ice  to  have  strength  values  over  twice  that  of  fresh  v/ater  ice 
(line  Dj-D),  provided  that  identical  testing  procedures  are  used.  This 
increase  in  strength  is  dependent  upon  the  amount  of  solid  salt  that  has 
precipitated.  Later  tests  have  shown  that  these  values  are  excessive. 

When  this  postulated  behavior  is  plotted  on  a  graph  that  shows  the  relative 
strength  (the  strength  of  sea  ice/the  basic  strength  of  sea  ice  o-0 )  vs  0^, 
it  appears  as  shown  in  Figure  21.  Assur  also  postulated  other  possibilities 
including  a  simple  curve  a(  versus  \Tv {^convex  to  the  center  of  coordinates. 
This  now  appears  to  be  more  justified. 


Figure  20.  Ring  tensile  strength  of  sea  ice  vs  square  root  of  the  relative 
brine  volume  vjj.  Each  point  represents  on  the  average  nine  tests  (Assur). 17 
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Figure  21.  Relative  ring  tensile  strength  of  sea  ice  (<rf/or0)  as  a  function  of 
temperature  and  salinity.  The  righthand  scale  indicates  a  strength  increase 
of  1/3  associated  with  the  precipitation  ofNaz  SO*  •  10H2  O.  Assur17  recom¬ 
mended  this  scale  be  used  for  perennial  ice  and  for  annual  ice  in  the  tem- 
•  perature  range  -8.  2  >  0  >  -22.  9C. 


Soon  after  the  publication  of  this  paper,  Langleben85  reported  on  the 
ring  tensile  strength  of  ice  fom  Shippegan,  New  Brunswick,  Canada.  He 
considered  only  ice  below  the  transition  zone  which  in  this  case  was  located 
approximately  23  cm  below  the  upper  ice  surface  because  of  the  formation 
of  an  appreciable  thickness  of  infiltrated  snow  ice.  No  pronounced  correla¬ 
tion  between  Of  and  either  or  \Tv  was  found.  He  did,  however,  note 

an  increase  in  Of  with  increasing  depth  in  the  ice  sheet  z  and  suggested 
that  this  effect  was  the  result  c  brine  between  the  crystals  being  more 
efficient  in  reducing  strength  than  brine  within  the  crystals  (the  micro¬ 
structure).  This  argument,  however,  seems  improbable  because  of  the 
extremely  small  volume  "sampled"  by  the  ring  tensile  test.  Even  in  the 
small  crystals  just  below  the  transition  layer,  the  average  grain  size  d  is 
four  times  the  size  of  the  largest  dimension  of  the  ring  tensile  test  volume. 
Therefore  in  most  cases  the  ring  tensile  test  would  not  "see"  the  grain 
boundarie  s. 
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In  an  attempt  to  resolve  a  num¬ 
ber  of  the  questions  posed  by  Assur's 
paper, 17  Weeks18*  ran  an  extensive 
series  of  tests  on  ice  tormed  by 
freezing  NaCi  solutions.  He  found 
a  simple  linear  relation  between  tr^. 
and  s/vjj (Fig.  22)  at  temperatures 
warmer  than  -21 . 2C  with  a  af  axis 
intercept  <r0  =  24.7  kg/cm2  and  a 
v0  value  of  231  %•.  The  value  of  cro 
was  slightly  lower  than  the  af  value 
determined  by  similar  runs  on  fresh 
water  ice  (29.  6  kg/cm2  at  -IOC). 

This  suggests  that  the  ratio  of  stress 
concentration  in  salt  ice  to  that  in 
fresh  water  ice  is  in  the  range  1.  1 
to  1. 5  and  defeats  the  previous  hy¬ 
pothesis  of  Anderson  and  Weeks4 
who  assumed  that  the  maximum 
strength  of  sea  ice  was  one-third  that 
of  fresh  ice  based  on  the  assumption 
of  stress  concentrators  in  sea  ice. 

Tests  on  NaCl  ice  colder  than  the 
eutectic  temperature  indicate  that 
the  strength  is  essentially  independent 
of  both  the  sample  temperature  and  the 
volume  of  NaCl*  2H2  O  in  the  sample.  When  the  strength  of  the  ice  matrix  is 
corrected  to  -  IOC,  the  strength  of  the  NaCl*  2H2  O  ice  is  26.5  kg/cm2  or  only 
slightly  more  than  <ro  .  There  is,  therefore,  no  indication  of  a  pronounced 
strengthening  effect  associated  with  the  precipitation  of  NaCl*  2HZ  O.  Assur's 
data  were  also  studied  relative  to  the  possible  1/3  increase  in  strength  asso¬ 
ciated  with  the  precipitation  of  Na2  S04  •  10H2  O.  This  reanalysis  suggests 
either  (1)  that  the  effect,  if  it  exists,  causes  a  strength  increase  of  20%  or 
less,  or  more  likely  (2)  that  there  is  no  strength  increase  associated  with 
Na2  S04  •  10H2  O.  Peyton121  has  again  recently  discussed  possible  Na2  SO4  • 

10H2  O  reinforcement  on  the  basis  of  compression  tests.  These  results  will 
be  discussed  iater. 

Several  experiments  were  made  to  check  the  suggestion  of  Butkovich14 
that  the  change  between  the  <r£  values  at  Hopedale  and  Thule  was  produced  by 
differences  in  the  past  thermal  history  of  the  sample.  It  was  found  that  o-£ 
for  samples  tested  at  temperatures  higher  than  the  eutectic  temperature  can 
be  considered  as  independent  of  the  details  of  the  thermal  history  immediately 
before  testing  (Fig.  23).  This  assumes,  of  course,  that  extensive  brine  drain¬ 
age  or  redistribution  has  not  occurred.  On  the  ether  hand,  the  strength  of  ice 
tested  at  temperatures  below  the  eutectic  point  is  strongly  dependent  upon 
whether  or  not  the  specimen  has  been  subjected  to  temperatures  above  the 
eutectic  point.  This  loss  of  strength  is  presumably  due  to  brine  drainage  or 
redistribution  in  the  small  near  surface  volume  of  ice  where  failure  is  initiated. 
It  may  also  result  from  microcracks  introduced  by  the  phase  changes.  There 
was  no  evidence  for  any  significant  metastable  lag  between  the  time  a  sample 
reaches  a  given  temperature  and  the  time  the  "equilibrium"  cr£  value  is  reached. 
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Figure  22.  Average  ring  tensile 
strength  vs  square  root  of  the  av¬ 
erage  brine  volume,  NaCl  ice 
(Weeks). 187 


34 


THE  MECHANICAL  PROPERTIES  OF  SEA  ICE 


30 


20 


10 


r~ 


[UU\_ 


1  1  1  I  1  1 


2  3  4 

TIME  (OATS) 


-J 

5  6 


Figure  23.  Study  of  the  effect  of  ther¬ 
mal  cycling  on  the  ring  tensile  strength 
of  NaCI  ice;  test  temperature  -20C, 
cycled  to  -IOC  (Weeks).  187 


Figure  24.  Variation  in  brine 
volume  v0  ,  producedina  con¬ 
stant  width  model  (d0  =  7xl0-3 
cm)  as  a  result  of  the  varia¬ 
tion  of  50  with  z  (Weeks  and 
Assur).  187 


In  a  further  analysis  of  this  NaCI  ice  data,  Weeks  and  Assur189  pointed 
out  that  a  systematic  increase  in  the  value  of  the  plate  spacing  a0  with 
increasing  depth  in  the  ice  sheet  z  should  result  in  a  systematic  decrease  in 
the  /alue  of  v0  .  Using  the  a0  =  f(z)  relation  observed  for  this  NaCI  ice  (Fig. 
12)  and  assuming  a  constant  width  model  with  d0  =  7  x  10'3  cm,  the  calculated 
change  in  v0  vs  z  is  shown  in  Figure  24.  They  then  examined  the  NaCI  ice 
ring  tensile  results  to  see  if  such  a  change  could  be  found.  Figure  25  shows 
their  results:  a  systematic  decrease  in  v0  as  determined  from  <rf  vs  \Tv~b 
plot 3  as  z  increases.  They  then  suggest  that  similar  results  should  occur  in 
natural  sea  ice.  They  explain  the  pronounced  increase  in  irf  with  increasing 
?.8S  as  the  result  of  the  deterioration  in  the  strength  characteristics  of  the 
upper  part  of  the  ice  sheet  during  the  formation  of  infiltrated  snow  ice. 
Further  theoretical  considerations  based  on  these  results  are  discussed  in 
Assur  and  Weeks30  and  were  outlined  earlier  in  this  paper. 

Graystone  and  Langleben57  reported  on  results  from  annual  ice  in  the 
Canadian  Arctic.  Their  tests  were  performed  immediately  after  the  sample 
was  removed  from  the  ice  sheet  and  the  resulting  equation  was 


rrf  =  29.0 


V 


yb 

o.  296 


(4.3) 


for  samples  whose  ice  temperature  was  higher  than  -8.2C.  This  was  ad¬ 
ditional  strong  experimental  evidence  in  favor  of  the  previously  assumed  sea 
ice  models  and  was  also  in  good  agreement  with  the  results  of  Weeks187  on 
NaCI  ice.  Note  that  the  <ro  value  is  almost  identical  with  the  ring  tensile 
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strength  of  pure  ice  as  determined  by 
Weeks.  The  ov  values  for  samples 
colder  than  -8.2C  and  therefore  con¬ 
taining  Na2  SO*  •  10H2  O  fall  on  this 
same  curve.  These  authors  obtained 
very  poor  results  when  samples  stored 
at  temperatures  warmer  than  -20.  5C 
were  used.  These  erratic  results  are 
presumably  produced  by  brine  drainage. 
Although  brine  drainage  corrections 
were  made,  no  significant  improvement 
in  the  quality  of  the  data  was  noted. 

Results  similar  to  those  of 
Graystone  and  Langleben  were  obtained 
by  Dykins44 


crf  =  30.  14 


1  - 1 

(1  Vb 

y  0. 585 

(4.4) 


Figure  25.  Decrease  in  the  brine 
volume,  v0  (%o),  with  increasing 
distance  below  the  surface  of  the 
ice  sheet,  z.  The  v0  values  were 
determined  from  plots  of  j-f  vs 
(Weeks  and  Assur).  ^ 


on  ice  cylinders  in  which  no  holes  had 
been  drilled.  Dykins  assumed  that 
K  =  6  (that  brine  pockets  were  acting 
as  infinitesimally  small  holes)  in 
making  these  calculations.  But  using 
the  empirical  K  of  5.2  determined  by 
Frankenstein50  we  obtain  a  (r0  value 
of  27.2  kg  /cm2  for  Dylans'  data. 

Langleben  and  Pounder88  expanded 
their  studies  to  two  year  old  and  polar 
(multi-year)  ice  at  Isachsen,  N.W.T. 
Because  of  low  air  temperatures  .  -10  to  -15C)  and  low  ice  salinities  (Sj<l%o), 
'Jv^was  always  less  than  0.  1  (vjj  <10%o).  Graystone  and  Langleben,  57  on  the 
other  hand,  hadv^  values  up  to  100%?.  Comparing  their  results  with  the  strength 
of  annual  ice  with  the  same  brine  volume  (eq  4.  3),  they  found  that  biennial  ice  is 
21%  stronger  than  annual  ice  and  6%  stronger  than  polar  ice.  These  differences 
are  difficult  to  explain,  Assur17  also  found  higher  strength  values  for  warm 
perennial  ice.  Hendrickson  and  Rowland's62  data  for  perennial  ice  are  also 
slightly  higher.  A  slight  but  apparently  non- significant  difference  was  found 
between  o-f  values  from  horizontal  and  vertical  cores:  <rf  (horizontal)  >  crf  (ver¬ 
tical  by  4%.  It  should  be  noted  that  the  horizontal  cores  were  not  systematically 
placed  in  the  'hard  fail"  orientation  in  the  manner  used  by  Butkovich36  who 
found  that  horizontal  cores  gave  consistently  higher  values  than  vertical  cores 
from  annual  ice.  They  observed  that  the  standard  deviation  of  the  <rf  values  for 
horizontal  cores  (7.6  kg/cm2  )  was  appreciably  greater  than  that  of  the  vertical 
cores  (4.7  kg/cm2  ). 

The  most  extensive  study  of  ring  tensile  results  currently  available  is  that 
of  Frankenstein50  who  reports  the  results  of  over  1400  individual  tests  (Fig.  26), 
Each  point  in  this  figure  represents  roughly  nine  tests.  The  least- squares 
equation  for  a  straight  line  fit  to  the  data  when  s/v^  <  0.400  is 


af  =  28.51 


■f 

"0. 
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(4.5) 
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Figure  26.  Ring  tensile  strength  vs  square  root  of  the  brine  volume,  sea 

ice  (Frankenstein). 50 


Equation  4.5  is  in  good  general  agreement  with  the  results  of  Weeks187  and 
Langleben  and  Pounder.88  When  >  0.400,  o-f  remains  constant  at 
approximately  6.7  kg/cmz  .  It  is,  however,  very  doubtful  that  this  is  the 
proper  relation  at  high  brine  volumes  inasmuch  as  the  equation  4.  1  used 
to  calculate  these  results  assumes  a  high  stress  concentration  near  the  hole. 
Plastic  relief  undoubtedly  takes  place  at  high  brine  volumes;  therefore  the 
true  relation  in  this  range  will  probably  lie  between  the  laboratory  results  of 
Weeks  and  the  field  results  of  Frankenstein.  A  possible  way  of  incorporating 
the  o-f  values  at  * Jv ^  >  0.  400  into  a  general  theory  was  discussed  earlier. 
Figure  26  does  not  show  an  appreciable  strength  difference  in  the  Na2  S04  • 
10H2  O  range. 

Although  there  is  no  obvious  reason  to  expect  Antarctic  sea  ice  to  be 
different,  it  is  interesting  to  investigate  this  point.  Figure  2762  shows  the 
arctic -ice  curve  of  Frankenstein,  an  Antarctic -ice  curve  of  Hendrickson  and 
Rowland  which  includes  results  from  one-and  two-year  ice  and  perennial  ice, 
and  a  plot  of  the  1964-65  data  from  McMurdo  Sound  and  Hallett  Station.  Al¬ 
though  small  differences  in  the  slope  remain  to  be  explained,  the  curves  are 
in  reasonable  agreement  indicating  a  a  0  value  of  roughly  29  kg/cm2  .  Un¬ 
fortunately  there  are  not  sufficient  Tf  values  with  */vjj  >  0.400  to  prove  that 
( j-j.  ~  6.  7  kg/cm2  at  high  brine  volumes.  The  data,  however,  give  no  reason 
to  suspect  that  such  is  not  the  case.  Other  tests  available  from  the 
Antarctic  are  those  of  Abele  and  Frankenstein.1  Their  results,  however, 
cover  only  a  very  limited  range  of  values  making  estimates  of  <ro  and  v0 
uncertain. 

Smirnov15*  has  reported  the  results  of  a  large  number  of  "roller"  tensile 
tests.  This  test  is  performed  by  compressing  10  x  10  x  20  cm  rectangular 
blocks  of  sea  ice  between  two  cylindrical  rollers,  and  appears  to  give  slightly 
higher  strengths  than  the  ring  tensile  test.  The  results  clearly  show  a  pro¬ 
nounced  variation  with  changes  in  the  salinities  of  the  samples.  It  appears 
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Figure  27 .  Ring  tensile  strength  of  Antarctic  sea  ice  compared 
with  Arctic  sea  ice,  1964-65  (Hendrickson  and  Rowland).*2 

quite  probable  that  this  test  will  give  a  tr[  vs  'dvj,  relation  similar  to  that 
established  for  ring  tensile  tests. 

Direct  tensile  tests 

Until  recently  there  have  been  no  detailed  investigations  using  the  con¬ 
ventional  tensile  strength  test.  Fortunately,  Dykins4’  has  just  published  a 
series  of  direct  tensile  tests  performed  on  laboratory-grown  sea  ice.  The 
reduced  cross  section  in  his  tests  had  an  area  of  ~  13  cm2  and  the  extension 
rate  was  1.2  ern/min.  The  ice  showed  plate  spacing  values  that  increased 
from  0.  35  mm  near  the  top  of  the  ice  sheet  to  0.  45  mm  near  the  bottom. 

A  very  pronounced  increase  in  grain  size  was  noted  over  this  same  vertical 
distance,  from  8  grains/cm2  at  the  top  to  0.  3  grains/cm2  near  the  bottom 
(an  increase  by  a  factor  of  25  in  41  cm  of  ice).  Because  his  strength  profiles 
show  no  pronounced  vertical  variation,  this  can  be  considered  as  evidence 
that  changes  in  grain  size  have  no  major  influence  on  strength.  There  was  a 
consistent  variation  with  '/v"|-)  however  (Fig.  28),  of  the  same  gen?ra!  trend 
as  shown  by  the  ling  tensile  tests.  Although  the  data  are  limited,  there  does 
not  appear  to  be  any  significant  change  in  the  strength  characteristics  as¬ 
sociated  with  the  precipitation  of  solid  salts.  The  orientation  of  the  specimen 
relative  to  the  direction  of  growth  of  t.te  sea  ice  has  a  pronounced  effect  on 
tensile  strength.  As  should  be  expected  in  this  test,  vertical  specimens  which 
are  forced  to  break  across  the  ice  plates  show  strengths  similar  to  horizontal 
specimens  in  the  "hard  fail"  orientation  but  higher  strengths  (2  to  3  times) 
than  horizontal  specimens  in  the  "easy  fail"  orientation.  As  a  result  of  these 
differences,  horizontal  specimens  show  an  appreciably  greater  scatter  than 
vertical  specimens.  Dykins  suggests  that  the  results  of  these  direct  tension 
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Figure  28.  Tensile  strength  of  sea  ice 
vs  square  root  of  the  brine  volume.  In 
the  stronger  samples,  tensionis  exerted 
parallel  to  the  growth  direction  and  in 
the  weaker  samples,  perpendicular. 
The  brackets  indicate  the  95%confidence 
limits  on  the  mean  (Dykins).4* 

Prominent  among  the  earlier  studie 
who  on  the  basis  of  least-squares  gave 

cr,  =  4.7  -  0.96  0  -  0. 31  02 


tests  can  be  compared  directly  with 
flexural  strength  results  for  sea  ice. 

If  further  tests  substantiate  this 
suggestion,  this  type  of  testing  pro¬ 
cedure  should  prove  to  be  very  use¬ 
ful. 

The  only  other  direct  tension 
results  are  those  of  Peyton. ,22»123 
Unfortunately  these  results  are 
currently  presented  in  a  form  which 
makes  them  difficult  to  use  or  com¬ 
pare  with  other  published  results. 

Flexural  Strength 

Small-beam  tests 

The  small-beam  flexural  strength 
test  was  among  the  tests  first  applied 
to  sea  ice  studies.®0’*5  Recently, 
however,  the  test  appears  to  have 
lost  "favor"  relative  to  the  ring  tensile 
test.  The  reasons  for  this  are  two¬ 
fold.  Ring  tensile  samples  are  easier 
to  prepare  and  also  give  better  correla¬ 
tions  with  such  parameters  as  brine 
volume  and  microstructure  than  do 
small  sample  flexural  tests. 

i  were  the  results  of  Arnol'd-Aliab'ev11’ 15 


(4.6) 


where  af  is  in  kg/cm2  ,  as  the  temperature  dependence  (-O.4>.0>,  -18.  5C)  for 
small  beam  flexural  results.  He  then  interpreted  the  deviations  from  this  curve 
in  terms  of  the  salinity  of  the  ice  Sj  (%o)  which  in  turn  he  empirically  related  to 
air  porosity  Va(cm3/kg) 

V  =  96.  -  S..  <4-7) 

a  l 


Combining  his  equations  we  obtain 

ff,  =  +  18.76  -  0.96  0  -  0.31  02  -  0.15  V  (4.8) 

f  a 

which  indicates  that  <j-f  increases  as  0  and  V  decrease.  It  is  interesting  to 
note  that  although  Arnol'd-Aiiab'ev  pioneered  techniques  for  the  determination 
of  the  air  content  of  sea  ice,  he  correctly  realized  that  a  considerable  per¬ 
centage  of  the  air  porosity  is  the  result  of  brine  drainage  prior  to  testing. 

The  first  set  of  recent  small  beam  results  accompanied  by  sufficient 
supplementary  information  for  an  analysis  was  published  by  Butkovich34’ 3* 
who  performed  tests  on  both  horizontally-,  and  vertically-,  cut  center-loaded 
simple  beams.  In  the  horizontal  beams,  the  direction  of  loading  was  normal 
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Figure  29.  Small  beam  flexural  strength  of  sea  ice  vs  brine  volume 

(Butkovich). 36 

tc  the  ice  sheet  while  in  the  vertical  beams  the  long  axes  of  the  crystals 
were  parallel  to  the  length  of  the  beam,  The  test  results  are  shown  in 
Figure  29.  There  is  a  large  scatter  and  vertical  beams  give  significantly 
higher  values  than  horizontal  beams.  No  appreciable  difference  is  indicated 
between  the  Hopedale  and  Thule  results.  This  is  interesting  inasmuch  as  a 
definite  difference  was  indicated  by  ring  tensile  tests.  Note  also  that  the 
small  beam  values  are  consistently  lower  than  ring  tensile  results  from  the 
same  ice  (8  kg/cm2  as  compared  with  18  kg/cm2  ).  The  pronounced  de¬ 
pendence  of  irf  on  Vfc  as  shown  in  ring  tensile  tests  is  also  absent.  Quantita¬ 
tively  similar  results  have  been  obtained  by  Petrov121  and  Smirnov, 153 
Tabata,!62  and  by  Abele  and  Frankenstein1  for  ice  from  the  Arctic  Ocean, 
the  Okhotsk  Sea,  and  from  McMurdo  Sound  in  the  Antarctic  respectively. 
Qualitatively  similar  results  were  also  reported  by  Pounder  and  Little125 
although  no  salinity  determinations  were  made  on  their  samples.  Russian 
observers145,  147  *  24  have  also  made  similar  measurements  in  the  vicinity  of 
Mirny  Station  in  the  Antarctic.  Unfortunately  the  only  supplementary  data 
they  provide  are  air  temperatures145, 147  and  salinities  and  densities.24 
They  do  show,  however,  that  during  the  testing  periods  the  maximum  ice 
strength  occurred  near  the  center  of  the  ice  sheet, 

Butkovich35  has  attempted  to  explain  the  observation  that  ring  tensile 
tests  give  larger  values  than  small  beam  and  in- situ  flexural  tests  on  the 
following  basis.  If  it  is  assumed73  that  sea  ice  contains  a  distribution  of 
flaws  and  that  each  of  these  flaws  can  withstand  stresses  up  to  a  certain 
size,  then  it  follows  that  the  probability  of  the  "weakest"  imperfection 
occurring  in  this  volume  is  proportional  to  the  volume  size.  Therefore  tests 
that  force  failure  on  a  small  volume  (ring  tensile)  should  presumably  yield 
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higher  strength  values.  As  Goetze53  has  pointed  out,  this  volume  effect  is 
very  pronounced  in  coal  and  glass  where  the  flaw  distribution  shows  a  wide 
scatter.  On  the  other  hand,  in  polycrvstalline  materials,  where  the  flaws 
are  related  to  the  microstructure,  the  variation  in  well-annealed  samples 
is  much  less.  A  reanalysis  by  Goetze  of  Jellinek's73  results  concerning 
the  effect  of  the  sample  volume  on  the  strength  of  fresh  water  ice  shows 
that  the  strength  appears  to  be  independent  of  volume  as  long  as  the  mini¬ 
mum  sample  dimension  is  larger  than  roughly  2  mm.  This  implies  that  in 
samples  with  dimensions  smaller  than  2  mm  there  is  some  "interference" 
with  the  microstructure  that  is  initiating  the  failure.  It  should  be  recalled 
that  earlier  arguments  (eq  3.  48)  suggested  that  in  sea  ice  the  size  of  this 
microstructure  is  in  the  range  0.  1  to  1.0  mm. 

Goetze  further  argues  that  the  difference  between  the  small  beam 
flexural  and  the  ring  tensile  values  is  caused  by  different  aspects  of  the 
microstructure  controlling  the  fracture  mechanism.  In  small  beams  the 
fracture  may  be  related  to  the  dimensions  of  the  crystals  as  compared  to 
the  dimensions  of  the  brine  pockets  in  ring  tersile  tests.  Some  justification 
for  this  suggestion  may  be  obtained  by  calculating  the  constant  in  eq  3.48 
from  Butkovich's37  results  on  the  flexural  strength  of  well-annealed,  uniform, 
isotropic  glacier-ice  from  the  Tuto  Tunnel.  When  Goetze  does  this  he  ob¬ 
tains  an  average  value  for  (<rf  d)  of  111.  Now  extending  this  argument  to  the 
typically  larger  grain  diameters  encountered  in  sea  ice  (1  to  2  cm)  he  obtains 
<r £  values  of  10.5  and  7. 5  kg/cm2  respectively.  These  values  are  in  the 
range  of  values  as  determined  by  small  beam  flexural  tests.  It  should  be 
quite  simple  to  verify  this  suggestion  of  a  grain  size  related  fracture  mech¬ 
anism.  All  that  is  required  is  a  series  of  small  beam  tests  on  ices  with  a 
constant  brine  volume  and  a  range  of  grain  sizes.  It  shoul  i  be  pointed  out 
that  Dykins'45  direct  tension  results  do  not  appear  to  support  the  grain  size 
argument.  If  the  fracture  mechanism  was  related  to  grain  size  in  Dykins' 
tests,  a  pronounced  vertical  variation  should  have  been  observed  in  the 
failure  strength.  This  was  not  the  case.  Because  of  the  large  cross  sectional 
area  of  his  tensile  samples  (~13  cm2  ),  the  test  volume  does  not  appear  to  be 
of  a  size  that  would  interfere  with  a  grain  size  related  fracture  mechanism. 
Another  reason  for  the  reduced  strength  of  larger  samples  is  the  presence  of 
large  holes  and  defects  such  as  brine  drainage  channels25  and  macro-cracks. 
The  probability  that  such  defects  would  occur  in  a  small  sample  is  quite  low. 

On  the  other  hand,  ring  tensile  tests  may  be  subjected  to  an  unknown 
plastic  reduction  in  stress  concentration  as  suspected  by  Sala.*37  Evidence 
for  this  suggestion  can  be  found  in  the  tests  of  Smirnov153  who  compared  a 
large  number  of  "roller"  tensile  values  with  flexural  results  determined  on 
similar  ice.  His  results  are  shown  in  Figure  30.  Theoretically  the  flexural 
strength  should  be  equal  to  the  tensile;  in  practice  this  is  rarely  observed. 
Nevertheless  it  might  be  expected  that  the  ratio  of  the  two  strengths  would  be 
constant.  Figure  30  shows  that  this  is  not  the  case.  For  low  salinities  sea 
ice  is  brittle  and  the  tensile  and  flexural  strength  is  roughly  equal.  For  high 
salinities  sea  ice  is  plastic  and  the  tensile  strength  is  considerably  higher 
than  the  flexural  strength.  We  suggest  that  this  is  caused  by  the  plastic 
relief  of  stress  concentrations  near  everpresent  small  holes  in  the  "roller" 
tensile  samples.  Smirnov's  data  therefore  indicate  that  plastic  relief  ol 
stress  concentrations  may  be  an  important  factor  in  certain  types  of  tests. 

The  small  beam  flexural  tost  has  also  been  extensively  utilized  by 
Tabata152  ’  158  to  study  the  relat’on  between  failure  strength  and  rate  of 
applied  stress.  These  results  will  be  discussed  later.  A  reliable  approach 
to  evaluate  the  strength  of  sea  ice  sheets  from  small  sample  tests  is  not 
available  as  yet. 
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Figure  30.  Average  values  of  "roller" 
tensile  strength  vs  small  beam  flexural 
strength  as  determined  on  similar  sea 
ice  samples  under  similar  conditions 
(Smirnov).  153 


In-situ  cantilevers 

The  in-situ  cantilever  or  key 
test  was  first  described  by  Neronov106 
in  conjunction  with  a  series  of  tests 
on  lake  ice.  Although  cantilever  re¬ 
sults  have  commonly  been  used  in 
estimating  erf  values  in  bearing  capac¬ 
ity  problems  via  relations  such  as 
eq  1.  12,  there  is  surprisingly  little 
documentation  of  the  results  of  such 
tests  in  the  literature.  The  reasons 
for  this  paucity  of  data  are  obvious 
to  anyone  who  has  ever  performed  or 
tried  to  analyze  the  results  of  canti¬ 
lever  tests.  The  test  is  quite  time- 
consuming  if  the  ice  is  thicker  than 
~50  cm.  Therefore  at  the  end  of  a 
testing  program  one  invariably 
finishes  with  an  insufficient  amount 
of  data  to  "make"  a  satisfactory 
technical  report.  In  addition  the 
analysis  of  cantilever  tests  is,  as  is 
the  analysis  of  any  test  that  considers 
the  ice  sheet  as  a  whole,  far  from 
unambiguous. 


The  only  extensive  study  utilizing  the  results  of  cantilever  testa  is  that 
of  Weeks  and  Anderson185  who  performed  208  tests  on  thin  (less  than  40-cm) 
sea  ice.  They  demonstrated  the  dependence  of  the  flexural  strength  on  the 
average  ice  salinity  and  temperature  but  did  not  plot  their  results  against 
brine  volume.  This  is  done  in  Figure  31.  Figure  31  also  contains  the  re¬ 
sults  of  Butkovich34  and  those  of  Brown32  (his  Figure  11)  that  have  not 
been  subjected  to  high  temperature  cycling.  The  equation  obtained  from 
Figure  31  is 


for  n/v>,  <  0.  33.  When  >  0.  33,  cf  is  constant  with  a  value  of 

roughly  2.0  kg/cmz  .  Although  additional  data  are  badly  needed,  these  re¬ 
sults  are  quite  encouraging  in  that  the  strength  variation  appears  similar  to 
that  found  by  ring  tensile  tests  although  strength  values  are  much  lower. 

Note  also  that  the  v<j  values  are  quite  similar:  202  %o  as  compared  with 
234  %o  . 

In  these  tests  little  difference  was  noted  between  "pull-ups"  and  "push¬ 
downs.  "  This  is  in  contrast  to  results  on  lake  ice  where  the  "pull-up" 
strength  is  significantly  higher  than  the  "push-down"  strength. 48*  49  This 
difference  is  presumably  the  result  of  the  presence  of  extensive  thermal  sur¬ 
face  cracks  on  lake  ice.  55 

It  is  interesting  to  note  the  pronounced  effect  that  a  warm  temperature 
cycle  had  on  Brown's32  cantilever  results.  This  is  in  striking  contrast  to 
the  experimental  results  of  Weeks187  where  no  change  was  noted  in  ring- 
tensile  specimens  cycled  and  tested  above  the  NaCl*  2H2  O  precipitation 
temperature. 
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Figure  31.  Flexural  strength  measured  by  in-situ 
cantilever  beam  tests  vs  square  root  of  the  brine 

volume. 


There  has  been  considerable  discussion  regarding  possible  stress  con¬ 
centration  at  the  sharp  corner  at  the  butt  of  the  beam  during  cantilever  tests. 
Comparisons  of  in-situ  cantilevers  with  in-situ  simple  beams  (Frankenstein, 
personal  communication)  suggest  that  stress  concentration  is  of  the  magni¬ 
tude  of  10%  or  less.  Similar  conclusions  have  been  proposed  on  the  basis 
of  limited  data  by  Sala.137 


Shear  Strength 

The  information  available  on  the  shear  strength  of  sea  ice  is  extremely 
limited.  In  addition  many  of  the  values  described  as  shear  strengths  are 
actually  the  result  of  mixed  mode  failures.  34,144  The  test  results  of 
Butkovich34  from  Hopedale,  Labrador,  were  obtained  by  encasing  the  core 
in  three  close  fitting  metal  cylinders  and  shearing  out  the  central  third  of 
the  core.  These  results,  when  plotted  against  v/v^,  show  a  large  scatter 
and  indicate  an  average  shear  strength  of  ~21  kg/cm2  for  n/v^  =  0.  19. 
There  is  no  obvious  correlation  between  o-f  and  v/v^  .  The  data  do  suggest 
a  slight  increase  in  <j-£  with  decreasing  test  temperature.  These  shear 
strengths  are  considerably  higher  than  the  tensile  strengths  reported  by 
Butkovich  on  vertical  cores.  The  reason  for  this  is  presumably  that  the 
shear  failures  are  forced  to  occur  across  the  columnar  grains  and  therefore 
must  break  across  the  ice  platelets. 

Serikov144  has  reported  the  results  of  11  shear  or  punch  tests  in  which 
failure  mode  is  questionable.  His  average  shear  stress  is  8.6  kg/cm2  with 
a  maximum  (11.2  kg/cm2  )  near  the  center  of  the  ice  sheet.  The  only 
supplementary  information  provided  is  that  the  air  temperature  varied  be¬ 
tween  -0.9  and  -4.  1C  during  the  tests.  Petrov121  has  also  reported  the  re¬ 
sults  of  a  large  number  of  shear  tests  on  arctic  pack  ice.  Unfortunately 
his  results  have  not,  at  present,  been  analyzed  in  terms  of  either  or  va. 
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The  best  set  of  shear  tests 
available  is  by  Paige  and  Lee114 
from  McMurdo  Sound.  Their  re¬ 
sults  (Fig.  32),  determined  on  7.  6 
cm  diam  cores,  were  obtained  by 
using  specially  designed  shearing 
heads  which  moved  at  20  cm /min 
during  the  application  of  the  failure 
load.  Note  that  their  shear  strengths 
are  consistently  lower  than  their 
ring  tensile  values  and  that  the 
shear  strength  decreases  as  s/v^ 
increases  in  the  same  general 
manner  as  in  the  ring  tensile  and 
in-situ  cantilever  results. 


Figure  32.  Shear  strength  (solid  cir¬ 
cles)  and  ring-tensile  strength  (open 
circles)  vs  square  root  of  the  brine 
volume  (Paige  and  Lee).114 


Compressive  Strength 


As  has  been  noted  by  several 
authors,  the  values  obtained  by 
crushing  tests  depend  to  a  con¬ 
siderable  extent  on  the  conditions 
of  deformation:  i.e.  the  dimen¬ 
sions  of  the  sample  and  the  stress  rate  or. 15,181  These  variations  are,  of 
course,  superimposed  on  the  variations  caused  by  differences  in  the  struc¬ 
ture  and  orientation,  the  specimen  temperature,  and  other  factors.  Unfor¬ 
tunately  there  have  been  consistent  differences  between  the  sample  geometry 
used  by  American  and  Canadian  as  compared  with  Russian  investigators. 

In  general  the  North  Americans  have  used  3  in.  (7.  6  cm)  diam  ice  cylinders 
with  a  3  to  1  ratio  of  length  to  width.  In  fresh  water  ice  as  the  length  to 
width  ratio  decreases,  the  strength  increases  due  to  different  failure  modes. 
Limited  data  bearing  on  this  problem  can  be  found  in  Butkovich.  35  Kudriavtsev82 
recommends  using  10-cm  cubes  while  Serikov146  has  used  4.  3-cm  cubes. 


Butkovich34-  36  found  that  the  compressive  strength  of  vertical  cores  of 
sea  ice  is  considerably  greater  than  that  of  horizontal  cores.  This  difference 
is  quite  reasonable  in  view  of  the  structure  of  sea  ice.  In  vertical  cores  fail¬ 
ure  is  forced  to  occur  across  the  ice  platelets,  while  in  horizontal  cores 
fracture  can  more  readily  take  advantage  of  the  planes  of  weakness  in  the  ice 
crystals.  Similar  differences  were  noted  by  Bruns  and  Deriugin198  who  also 
found  that  the  ratio  of  the  vertical  to  the  horizontal  compressive  strengths 
decreased  and  approached  unity  for  deteriorated  summer  ice.  For  verti¬ 
cal  specimens,  Butkovich  found  median  strength  values  ranging  from  78 
kg/cm2  at  -5C  to  roughly  120  kg/cm2  at  -16C.  Average  values  on  horizon¬ 
tal  cores  in  the  same  temperature  range  vary  from  21  to  43  kg/ cm2  . 
Serikov146  obtained  values  between  10  and  60  kg/cm2  at  -5C  for  4.  3-cm 
cubus  collected  at  Mirny. 


The  most  complete  set  of  compressive  strength  observations  available 
is  reported  by  Peyton123  who  ran  tests  on  a  large  number  of  samples  of 
different  sea-ice  petrographic  types  at  various  orientations  and  stress  rates. 
He  found  a  marked  dependence  of  compressive  strength  on  sample  orientation 
as  shown  in  Figure  33.  In  this  figure  the  upper  curve  represents  compression 
results  while  the  lower  curve  represents  direct  tension  tests.  T’  j  ice  used 
in  these  tests  was  bottom  ice  with  a  grain  size  larger  than  the  dimensions  of 
the  specimens.  Therefore  each  sample  can  be  considered  a  "single"  crystal. 
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Figure  33.  Average  failure  strength  in  com¬ 
pression  (solid  circles)  and  in  direct  tension 
(open  circled)  vs  sample  orientation:  bottom 
ice,  -IOC  (Peyton).121  For  orientation  nota¬ 
tion,  see  text. 


In  this  ice  type  c-axis  orientations  in  the  ice  sheet  are  always  close  to 
horizontal.  The  loading  angle  notation  is  as  follows:  the  first  number  gives 
the  angle  between  the  axis  of  the  test  sample  and  the  vertical,  while  the 
second  number  gives  the  angle  between  the  sample  and  the  c-axis  of  the  ice 
crystal.  Note  that  the  ratio  of  the  strength  obtained  from  vertical  cores  to 
that  obtained  from  horizontal  cores  is  about  3:1,  in  agreement  with  the  re¬ 
sults  of  Butkovich.  36 

Peyton  noted  a  relation  between  the  compressive  strength  and  the  stresr 
rate  o’,  the  brine  volume  of  the  ice  sample,  and  the  position  of  the  specimen 
in  the  ice  sheet.  These  results  were  obtained  by  multiple  regression  tech¬ 
niques,  but  unfortunately  Peyton  does  not  give  all  the  regression  constants. 
Therefore  the  results  are  at  the  present*  only  of  use  in  a  qualitative  manner. 

Inasmuch  as  the  failure  strength  o’f  is  a  function  of  o’,  Peyton  has 
attempted  to  "remove”  this  variation  as  follows.  If  the  regression  equation 
is  of  the  general  form 

07  =  ao-bX,yi  X2y*  - Xyn  (4.1 


a  new  parameter  or  can  be  defined  as 


X  Yn. 
n 


This  parameter  is,  of  course,  related  to  the  original  o’f  value  but  unfortu¬ 
nately  it  is  currently  impossible  with  the  information  provided  by  Peyton  to 


greatly  expanded  version  of  Peyton's  report121  is  currently  (1967)  being 
published  by  the  Arctic  Environmental  Engineering  Laboratory,  University 
of  Alaska. 
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Figure  34.  <rp  from  compression  tests  vs 
square  root  of  the  brine  volume  (Peyton).123 
For  the  difference  between  the  solid  and  dashed 
lines  see  discussion  in  text. 


interpret  an  values  in  terms  of  actual  measured  compressive  strengths.  It 
would  have  been  more  useful  to  simply  correct  all  strength  values  to  a 
common  value  of  o'.  We  can,  however,  consider  <tr  as  a  strength  index.  If 
o-r  is  plotted  against  \/v^( Fig.  34)  and  a  least  squares  straight  line  is  fitted 
through  those  data  points  where  the  ice  was  at  a  higher  temperature  than 
-8.7C,  the  dashed  line  results.  Peyton  interprets  the  fact  that  ice  strengths 
determined  at  temperatures  below  -8.7C  lie  above  this  line  to  indicate  that 
Na2  S04  •  10  H2  O  precipitation  is  effective  in  strengthening  the  ice.  This 
may  well  be  the  ease.  There  is,  however,  another  possible  interpretation 
of  Figure  34.  Remembering  that  in  other  types  of  tests  (Fig.  26,  31,  32) 
erf  decreases  as  a  linear  function  of  \l vj, at  low  brine  volumes,  and  becomes 
approximately  constant  at  high  brine  volumes,  we  may,  as  a  hypothesis, 
assume  that  the  five  data  points  at  n/v ^  >  0.6  are  associated  with  this 
latter  type  of  behavior.  If  this  is  the  case  the  data  may  be  approximated  by 
the  solid  line  shown  in  Figure  34  and  it  may  not  prove  necessary  to  invoke 
Na2  S04  •  10  H2  O  reinforcement.  Examination  of  Peyton's  results  as  a  whole 
suggests  that  they  would  be  well  fitted  by  a  simple  curvilinear  relation  of  the 
general  form  of  eq  3.40.  Regardless  of  the  salt  reinforcement  question, 

Figure  34  is  quite  important  in  that  it  shows  that  compressive  strengths 
can  also  be  analyzed  using  equations  3.35  and  3,40. 

Elastic  Modulus 

Of  all  the  mechanical  properties  of  ice,  the  elastic  or  Young's  modulus, 

E,  appears  to  have  been  the  most  extensively  studied.  Experimental  results  for 
pure  ice  were  available  as  early  as  the  1820's  and  are  summarized  in  tabular 
form  by  Weinberg  et  al. 1,3  and  Voitkovskii. 181  There  is  a  very  considerable 
scatter  in  the  E  values:  3  x  1G9  to  l  x  1011  dynes/cm2  for  static  measurements 
and  6  x  1010  to  1  x  1011  dynes/cm2  for  dynamic  measurements.  Dynamic 
measurements  are  invariably  larger  than  the  static,  show  a  smaller  scatter, 
and  are  more  readily  reproducible.  The  reason  for  these  differences  can 
be  explained  as  follows:181  in  static  tests  which  measure  strain  after  the 
application  of  a  load,  the  deformation  can  be  measured  only  after  a  finite 
time  interval.  Because  of  the  viscoelastic  behavior  of  ice,  this  interval 
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is  invariably  sufficient  to  permit  viscous  as  well  as  the  desired  elastic  de¬ 
formation.  Therefore,  a  static  E  determination  which  is  specified  by  the 
measurement  of  total  deformation  does  not  characterize  the  resistance  of 
ice  to  an  instantaneous  elastic  deformation,  but  also  a  viscous  tin.e  de¬ 
pendent  deformation  after  a  specific  time  interval.  The  greater  the  stress, 
the  more  important  the  viscous  and  creep  components  become',  the  elastic 
component  is  directly  proportional  to  the  stress  while  the  creep  rate  is 
proportional  to  the  stress  to  some  power  (~2  to  4  at  high  stress  levels). 
Therefore,  the  value  of  E  calculated  on  the  basis  of  strain  measurements 
decreases  as  stress  increases.  That  this  holds  for  fresh  ice  has  been 
shown  by  the  results  of  Pir.egin181  on  river  ice.  It  is  to  be  expected  that 
these  arguments  are  even  more  applicable  to  sea  ice. 

The  value  of  E  determined  by  dynamic  methods  is,  on  the  other  hand, 
based  on  the  determination  of  either  the  rate  of  propagation  of  vibrations  in 
the  ice  or  the  natural  (resonant)  frequencies  of  different  types  of  vibrations. 
Here  the  displacements  are  extremely  small  and,  for  many  purposes, 
anelastic  effects  can  be  neglected.  The  results  of  dynamic  measurements 
of  E  are,  therefore,  reasonably  reproducible  and  are  usually  preferred  over 
static  measurements  of  E.  A  review  of  the  different  methods  of  determining 
E  for  sea  ice  is  given  by  Anderson.  3 

Dynamic  measurements 

In-situ  seismic.  Two  main  approaches  have  been  used  to  obtain  dy¬ 
namic  determinations  of  E  for  natural  sea  ice.  The  first  utilizes  seismic 
field  techniques  to  measure  the  velocity  of  propagation  of  P,  S,  surface  and 
fie  rural  waves  in  the  ice  sheet.  If  the  P  and  S  wave  velocities  and  the 
density  are  known,  all  the  other  elastic  constants  can  be  calculated  provided 
that  the  ice  sheet  is  assumed  to  be  homogeneous  and  isotropic.  These 
simplifying  assumptions  are  reasonably  adequate  for  lake  ice,  but  are  far 
from  correct  for  sea  ice.  A  general  theory  for  treating  elastic  wave  prop¬ 
agation  in  a  medium  consisting  of  n  layers  of  a  transversely  isotropic  ma¬ 
terial,  which  well  approximates  the  real  properties  of  sea  ice,  has  been 
developed  by  Ander son. 8’7'  8  Unfortunately,  published  results  of  the  de¬ 
tailed  application  of  these  theoretical  developments  to  seismic  measure¬ 
ments  on  sea  ice  are  still  not  available.  Another  difficulty  in  analyzing 
seismic  results  lies  in  relating  the  results  to  the  "state"  of  the  sea  ice. 

For  flexural  waves,  which  are  controlled  by  the  "average"  properties  of 
the  ice  sheet,  the  salinity  and  temperature  profiles  can  be  used  to  calculate 
an  average  brine  volume  or,  preferably,  an  effective  flexural  rigidity. 

Average  ice  sheet  properties  may  also  be  used  with  some  confidence  during 
the  summer  months  when  the  ice  sheet  is  fairly  close  to  being  both  isothermal 
and  isosaline.  However,  during  the  sea  ice  growth  season  it  i3  necessary  to 
use  the  average  properties  of  the  cold,  higher  velocity  upper  portion  of  the 
ice  sheet  to  analyze  the  velocity  variations  of  P  and  S  waves.  These  waves 
do  not  travel  directly  along  the  ice  surface  but  penetrate  to  a  depth  which 
depends  on  the  physical  properties  of  the  upper  portion  of  the  ice  sheet. 
Unfortunately,  the  exact  "dimensions"  of  this  high  velocity  channel  are 
somewhat  ar  bitrary.  3'31,  32  This  introduces  an  additional  degree  of  un¬ 
certainty  into  the  analysis  of  field  results. 

A  representative  sampling  of  the  results  of  seismic  field  determinations 
of  the  elastic  modulus  of  sea  ice  is  given  in  Table  I.  An  examination  of  this 
table  indicates  that  E  varies  from  1.7  to  5.7  x  1010  dynes/cm2  when  measured 
from  flexural  waves  (Ef)  and  1.7  to  9.  1  x  1010  dynes/r.m2  when  determined 
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Table  1.  Representative  immu  field  determinations  of  the  rla3»ic  psrsmrtrr*  of  »rs  »ce:  F^  and  Ej, 
from  body  wive  and  atr  coupled  wave  measurements.  G.  shear  modulus,  u,  Poisson's  ratio,  and  h,  ice 


Ice  type  and  location 


Fast  ice.  North  c  Baltic 
Sea’ 


(10>J  dy/  cm*  )  ( IQ10  J 


Fast  ;ce,  Barter  Island, 
Alaska 

Spring 

-  1 80 

s.  s 

Fast  ice,  Hopedale, 
Labrador 

Winter 

... 

9.  1 

“ 

Fast  ice,  Hopedale. 
Labrador 

Spring 

-- 

3.  3  -  4.  H 

' 

Pack  ice,  Arctic  Ocean 

-- 

... 

s.  o 

- 

Pack  ice.  Bering  Sea 

Spring 

-  00 

s.  6 

- 

Pack  ice.  1  year.  Arctic 
Ocean.  NP-4 

Spring 

100 

6.  0 

* 

Pack  ice.  perennial, 

Arctic  Ocean,  NP-4 

Spring 

"-*400 

6.  H 

- 

Fast  ice.  Thule  Greenland 

Fall 

IS 

s.  s 

- 

Fast  ue,  Thule.  Greenland 

Fall 

S2 

-- 

c, 

Fast  ice.  Thule.  Greenland 

W inter 

112 

-- 

s. 

Fast  ice,  Thule,  Greenland 

Winter 

122 

7.  3 

- 

Fast  ice,  Thule,  Greenland 

Spring 

127 

-- 

1 . 

Fast  ice.  Thule,  Greenland 

Spring 

lt*2 

3.  1 

Fast  ice.  Cape  Sc  turn  ot. 
USSR 

Spring 

1  10-210 

5.  7 

Fast  ice.  Hokkaido.  Japan 

W  i  nt  e  r 

14-26 

2.9 

Pack  ice.  Arctic  Ocean, 
Station  Alpha 

Fall 

S.  « 

Pa<  k  ue.  Antic  Ocean, 
Station  Alpha 

Winter 

-who 

7.8 

Pack  ice,  Arctic  Ocean, 
Station  Alpha 

Summer  -^ilO 

7.  2 

FaM  ice.  Wales.  Alaska 

Winter 

-- 

4.  6 

Fast  ice,  Wales,  Alaska 

Spring 

-- 

3.  1 

Fa s t  ice.  Thule.  Greenland 

W inter 

1  is 

7.  7 

Fast  ice,  Hokkaido,  Japan 

Winter 

21 

1 .  7 

Fast  ice,  2  year, 

Isachsen,  N. W  .  T. 

Spring 

170 

7.  2 

from  in- situ  body 

wave 

velocities 

(Eb>- 

Eb  is  i 

/cm*  1  ( 1 0'°  dy/cmf  1 


rUltll  m  dull 
r  thickness. 

Source  and 
reference  nj, 
78 


-- 

Pome 

-- 

- 

Pome 

-- 

.  29 

.  17 

2.  2 

-- 

30 

2.  S 

-- 

30 

> 

3 

1.9 

.  37 

3 

3 

3 

94 

1.  1 

.  18 

68 

2.0 

.  36 

66 

3.  0 

.29 

66 

2.0 

.  38 

66 

1.7 

.  31 

32 

1.  1 

.  2  S 

52 

2.  8 

.  38 

32 

.  6 

.  17 

69 

2.  9 

.  2  S 

129 

determined  on  similar  ice.  This  is  reasonable  inasmuch  as  Ef  is  representa¬ 
tive  of  the  whole  thickness  of  the  ice  sheet  while  the  value  of  El  is  controlled 
by  the  properties  of  the  high  velocity  channel  in  the  usually  colder  and  stronger 
upper  portion  of  the  ice  sheet.  The  values  of  E  show  a  systematic  variation 
throughout  the  year  from  low  values  for  young,  warm,  high-salinity  new  ice 
to  high  values  for  cold  winter-ice.  A  decrease  occurs  again  in  the  spring 
when  the  ice  temperature  rises.  This  decrease  is  clearly  shown  by  Listov's 
results118  on  ice  from  the  Laptev  Sea.  As  the  temperature  of  this  ice  in¬ 
creased  in  the  spring,  the  ice  thickness  decreased  from  2  m  to  1 . 4  m  and  the 
values  of  E  showed  a  pronounced  decrease  (Fig.  35).  One  of  the  major  con¬ 
tributing  factors  to  this  variation  in  E  is  the  annual  variation  in  V  ,  the 
velocity  of  the  longitudinal-plate  wave.  The  value  of  Vp  has  been  aetermined 
on  thick  (~3  m)  pack  ice  throughout  a  complete  year  by  Hunkins6*  and  has 
been  found  to  vary  from  a  low  of  2300  m/s  in  the  su.nmer  to  a  high  of  3200 
m/s  in  the  winter  (Fig.  36).  Similar  but  less  detailed  res  i It s  oi  VD  vaiues 
and  flexural  wave  dispersion  curves  have  been  obtained  on  annual  sea  ice.-’ 
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Figure  35.  Decrease  in  elastic  Figure  36.  Seasonal  variation  of  longitu- 
modulus  E  with  increasing  ice  'em-  dinal  plate  velocity  (Hunkins).66 

perature  as  determined  by  in-situ 
seismic  measurements;  Laptev  Sea, 
spring  (data  of  Listov).113 


These  changes  in  Vp  and  in  E  are  clearly  related  to  changes  in  the  tempera¬ 
ture  of  the  ice. 29,30  ’  w  Figure  37  shows  a  plot  of  E  vs  temperature  (mean  ice 
temperature)  for  annual  fast  ice  with  a  thickness  of  1.  3  to  2.  1  m  from  Cape 
Schmidt,  Siberia.94  The  marked  increase  in  E  with  decreasing  temperature 
is  presumably  related  to  the  brine  volume,  but  because  detailed  salinity 
profiles  are  not  given,  the  brine  volume  is  difficult  to  estimate.  Fortunately 
the  relation  between  seismic  determinations  of  E  and  the  brine  volume  vl 
in  the  ice  at  the  time  of  the  determination  have  been  studied  by  Anderson3^ 
at  Thule,  Greenland,  (Fig.  38)  and  by  Brown  and  Howick31  on  the  Bering  Sea. 
Anderson  assumed  that  the  high  velocity  "channel"  extended  to  the  depth  where 
a  rapid  increase  in  V)j  was  noted  in  the  v^,  profile.  Brown32  (his  Fig.  9)  has 
made  a  similar  analysis  of  field  data  collected  at  Wales,  Alaska,  and  Thule, 
Greenland.  In  specifying  the  properties  of  the  high  velocity  "layer,  "  the  brine 
content  of  the  coldest  20-cm  layer  was  used.  Considering  these  difficulties 
the  results  of  both  Anderson  and  Brown  are  in  remarkably  gcod  agreement 
indicating  a  "pure"  ice  intercept  of  roughly  9  to  10  x  1010  dynes/cm2  for  Vb  =  0 
and  a  rapid  dec  'ease  in  E  as  increases.  An  interesting  attempt  to  develop 
a  theory  relating  sea  ice  strength  to  its  elasticity  via  a  perforated  plate  model 
has  been  made  by  Brown.  32 

Small  specimen.  The  main  alternate  approach  to  the  determination  of  the 
elastic  modulus  of  sea  ice  utilizes  small  samples  which  presumably  are 
reasonably  isothermal  and  isosaline.  By  using  either  cold  rooms  or  "heated" 
rooms  at  a  field  site,  a  much  wider  range  of  temperatures  and,  therefore, 
brine  volumes  can  be  studied  than  is  possible  when  the  ice  sheet  is  treated  as 
a  whole.  Considerable  care  must,  how  ”er,  be  taken  to  prevent  brir.e  drain¬ 
age.  The  natural  resonance  frequenci*  s  or  the  velocities  of  propagation  of 
sonic  waves  through  the  specimens  are  determined  and  the  elastic  constants 
are  calculated  from  the  results.  For  the  purposes  of  this  particular  review, 
the  small-specimen  results  are  more  interesting  than  the  field  seismic  fe- 
sults  and  will  be  discussed  in  more  detail. 
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Figure  37.  Elastic  modulus 
E  and  shear  modulus  G  as 
determined  by  iri-situ  seismic 
measurements  vs  mean  ice 
temperature,  annual  fast  ice, 
spring,  Cape  Schmidt,  Siberia 
(Lin'kov).94 


Figure  38.  Elastic  modulus  E,  as  de¬ 
termined  by  seismic  measurements, 
vs  brine  volume  (Anderson).  3 


The  first  small- specimen  results  were  reported  by  Tabata,  1581  1W.I62 
who  obtained  the  natural  resonance  frequency  of  sea-ice  bars  with  dimensions 
of  roughly  2x3.5x35  cm.  Even  though  the  initial  ice  blocks  were  obtained 
from  young  high  salinity  sea  ice,  the  long  storage  periods  before  actual  testing 
permitted  extensive  brine  drainage  causing  sample  salinities  to  be  in  the  range 
0.2  to  3.2%o.  Therefore,  attempts  to  correlate  variations  in  E  with  changes 
in  were  unsuccessful.  Inasmuch  as  most  of  the  initial  brine  had  been  re¬ 
placed  by  air,  a  plot  of  E  vs  v  =  1*5  +  va  was  also  prepared  (Fig.  39)  where 
va  was  measured  directly.  This  figure  suggests  a  strong  linear  dependence 
of  E  on  v  of  the  general  form 

#-  =  l  -  7*  (4.12) 

Lo  vo 


where  E0  is  the  sonic  modulus  of  fresh  ice197  and  v0  is  the  v  axis  intercept 
(vo  ~180  -  220%o).  Infiltr<  >ed  snow  ice  appears  to  have  slightly  higher  values 
of  E  for  a  given  value  of  v  than  does  normal  sea  ice.  This  is  to  be  expected 
inasmuch  as  the  ordered  brine  arrays  in  normal  sea  ice  crystals  should  cause 
E  to  decrease  more  rapidly  with  increases  in  v ^  than  would  occur  with  a 
random  array  of  brine  pockets. 

When  the  variation  in  E  with  specimen  temperature  was  studied,  Tabata 
found  a  slight  increase  with  decreasing  temperature  at  temperatures  lower 
than  -10C.  There  was  also  a  pronounced  decrease  in  E  at  temperatures 
higher  than  -10C.  Because  this  decrease  has  not  been  observed  in  fresh  ict,197 
it  undoubtedly  is  the  result  of  the  rapid  increase  in  the  volume  of  brine  in  the 
ice  at  temperatures  near  the  melting  point.  There  are  slight  changes  in  slope 
in  some  of  the  E  vs  temperature  plots  at  temperatures  of  roughly  -8  and 
-23C.161  It  is  interesting  that  while  these  changes  are  noticeable  in  most  of 
the  infiltrated  snow  ice  samples,  they  can  only  be  noted  in  one  of  the  eight 
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Figure  39.  Elastic  modulus,  E,  vs  volume  of 
air  plus  brine,  v;  young  sea  ice  (Tabata).159 
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sea  ice  curves.1  These  slope  changes  are  probably  the  result  of  both  rapid 
changes  in  the  brine  volume  and  the  presence  of  solid  salts.  There  has 
not  been,  at  present,  any  attempt  to  separate  these  two  effects. 

The  first  published  results  of  small  specimen  determinations  of  E  on 
sea  ice  that  had  not  undergone  significant  brine  drainage  were  by  Pounder 
and  Stalinsky. 127  They  measured  transit  times  for  pulses  of  supersonic 
frequence  through  sections  of  sea  ice  cores  from  Barrow  Strait,  west  of 
Cornwallis  Island,  N.  W.  T.  Although  they  observed  an  appreciable  scatter 
in  the  data,  there  clearly  appeared  to  be  a  linear  relation  between  E  and 
ice  salinity  at  a  constant  temperature.  This  implies  a  linear  relation 
between  E  and  v^. 

These  studies  were  continued  by  Langleben88  and  Langleben  and 
Pounder87  on  sea  ice  from  several  locations  in  the  Canadian  Arctic  and 
Greenland.  Their  results  (averages  from  over  300  measurements)  for  cold 
arctic  sea  ice  are  shown  in  Figure  40  and  indicate  a  linear  relation  between 
E  and  vjj.  The  scatter  is  small.  The  results  of  tests  on  warm  polar  ice 
also  appear  to  be  linear  on  a  E  vs  vj,  plot  but  with  a  lower  intercept  (E0  ~ 
8.8  x  1010  dy/cm*  )  and  a  smaller  slope.  However,  as  pointed  out  by 
Langleben  and  Pounder,87  not  as  much  confidence  should  be  placed  on  these 
values  because  of  brine  drainage  and  rapid  changes  in  the  phase  diagram  for 
warm  sea  ice. 

Further  measurements129  on  Z  year  old  and  multiyear  pack  ice  (Sj  ~  1  %o) 
at  Jsacnsen,  N.  W.  T.  ,  have  shown  that  both  give  E  values  roughly  4  to  5% 
lower  than  values  determined  under  similar  conditons  from  winter  ice  with 
an  identical  brine  volume.  An  examination  of  the  data  also  gives  the  distinct 
impression  that  the  dependence  of  E  on  vu  has  decreased.  However,  as 
Pounder  and  Langleben  have  pointed  out  this  impression  may  be  caused  bv 
the  limited  range  of  values  that  were  available  during  this  test  series 
There  were  no  differences  noted  between  the  values  obtained  on  the  Z  year 
and  on  the  multiyear  ice.  This  suggests  that  the  decrease  in  E  is  in  some 
way  associated  with  effects  that  occur  during  the  first  summer  and  is  appar¬ 
ently  not  intensified  by  cycling  through  further  melt  periods.  It  is  also 
interesting  to  note88  that  an  appreciable  difference  was  found  between  the 
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Figure  40.  Elastic  modulus  of  cold,  arctic 
sea  ice  vs  brine  volume:  small  specimen 
tests  (Langleben  and  Pounder).87 


ring-tensile  strengths  of  these  two  ice  "types.  "  A  similar  linear  relation 
between  E  and  Vv^  has  been  obtained  by  Abele  and  Frankenstein  on  annual  sea 
ice  from  the  McMurdo  Sound  region,  Antarctica.  Serikov146  has  also  meas¬ 
ured  similar  small  sample  values  of  E  (3.9  to  5.  0  x  1010  dy/cm2  )  on  1  year 
spring  ice  in  the  Antarctic  as  has  been  reported  using  standard  field  seismic 
methods  on  similar  ice  in  the  Arctic  (Table  I). 

Static  measurements 

There  have  been  only  a  very  limited  number  of  determinations  of  E  by 
static  methods.  Butkovich34  determined  values  from  the  slopes  of  the 
straight  line  portion  of  stress-strain  curves  in  compression.  The  values  of 
E  ranged  from  0.  4  to  1.0  x  10*°  dynes/cm2  and  mean  values  systematically 
increased  (0.76  to  0.83  x  1010  dynes/cm2  )  as  the  specimen  temperature  de¬ 
creased  (-4  to  -18C).  His  values,  however,  were  affected  by  deformation  in 
the  loading  apparatus. 

The  same  experimenter36  later  determined  E  on  simple  small  beams. 
Although  there  is  considerable  scatter  in  the  results,  a  linear  decrease  in 
E  with  increasing  vb  is  indicated  as  in  the  dynamic  tests.  In  contrast  to 
the  dynamic  results, 129  tests  run  on  vertical  beams  gave  appreciably  higher 
results  thc'n  tests  on  horizontal  beams.  There  also  appears  to  be  a  decrease 
in  E  with  time  after  removal  of  the  specimens  from  the  ice  sheet.  E  values 
determined  b/  similar  tests  have  also  been  reported  by  Tabata.168  The  most 
important  aspect  of  his  results  deals  with  the  dependence  of  E  on  o-  (see  p. 

56).  As  was  expected  E  values  determined  by  static  methods  were  consistently 
lower  than  those  determined  by  dynamic  techniques. 


Shear  Modulus  and  Poisson's  Ratio 


Accurate  determinations  of  the  velocity  of  propagation  of  shear  waves 
in  sea  ice  are  much  rarer  than  determinations  on  compressional  wave  veloc¬ 
ities.  In  general  it  has  been  found  that  in  ice  impulsive  shear  waves  are  not 
generated  by  standard  non-directional  explosive  sources  and  that  second 
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arrivals  cannot  be  picked  with  the  accuracy  of  the  firs!  compressional 
arrival.  Therefore,  either  mechanical  devices  or  shaoed  charges  must 

be  used  to  generate  the  shear  waves.  Hunkinsw  also  found  that  it  was  diffi¬ 
cult  to  obtain  good  shear  wave  arrivals  at  distances  of  more  than  250  m  from 
the  wave  source. 


An  examination  of  the  shear  modulus  (G)  values  presented  in  Table  I 
and  Figure  37  shows  that,  on  the  average  G  is  38%  of  the  value  of  E  and 
varies  in  a  similar  manner  with  the  "state"  of  the  sea  ice.  The  reason  for 
the  similarity  in  the  variation  of  E  and  G  can  be  seen  by  examining  eq  4.  13. 


G 


E 

2(1  +  v)' 


(4.13) 


Although  Poisson's  ratio  p  varies  slightly,  its  variation  is  small  compared 
to  the  variation  in  E.  Therefore,  in  general  G  follows  E. 

There  is,  at  present,  insufficient  information  on  the  variation  of  p  with 
the  state  of  the  sea  ice.  In  most  cases  (Table  I),  it  appears  that  p  is  high 
for  young  saline  ice  and  decreases  as  the  ice  becomes  colder  and  the  brine 
volume  decreases.  Support  for  this  conclusion  is  found  in  the  data  of  Lin’kov94 
who  calculated  p  from  the  general  relation 


H  = 


(4.  14) 


wh'  ;c  Vs  and  Vp  are  the  velocities  of  the  shear  and  compressional  plate 
waves  respectively.  We  express  these  results  in  the  form 


Figure  41,  (p  -  1/3)  vs  mean  ice 
temperature,  based  on  Lin'kov's 
in-situ  seismic  determinations, 
Cape  Schmidt,  Siberia.94 


where  p0  is  Poisson's  ratio  extrapo¬ 
lated  to  0°C,  pj  is  the  final  Poisson's 
ratio  at  very  low  temperatures  (taken 
to  be  equal  to  that  of  fresh  ice  0.  33),u 
and  0O  is  a  characteristic  temperature. 
The  results  are  shown  in  Figure  41. 

The  final  least-squares  equation  is 


p  =  0.333 

+  6.  105  xlO'2  expf-  -^-2—).  (4.  16) 
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A  decrease  in  p.  with  decreasing  average  ice  temperature  is  clearly  indicated. 
Langleben  and  Pounder87  have,  on  the  other  hand,  reported  no  systematic 
variations  in  p  as  measured  by  small- sample  resonance  methods  with  changes 
in  salinity,  temperature  or  ice  type.  They  recommend  that  an  average  value 
of  0.  295  be  used  for  sea  ice.  Further  experimentation  will  be  necessary  to 
eliminate  the  lack  of  agreement  in  these  results.  A  relation  between  p  and 
brine  volume  or  overall  porosity  v  should  be  explored. 


Time -Dependent  Effects 


Stress  rate 

Tensile  strength.  As  was  pointed  out  earlier,  in  tests  concerned  with  the 
elastic  behavior  of  sea  ice,  if  the  rate  of  stress  application  £r  is  greater  than 
0.  5  kg/cm2  -  sec,  it  has  been  assumed  that  the  failure  strength  is  independent 
of  or.  Actually  most  of  the  tests  discussed  in  this  section  were  performed  at 
a  constant  rate  of  strain  j.  The  stress  rate  or  is  constant  only  if  Hooke's  law 
holds.  The  assumption  of  the  independence  of  erf  and  ir  was  based  on  the  re¬ 
sults  of  Jellinek73  who  showed  that  the  strength  of  snow  cylinders  that  failed 
in  tension  behaved  in  this  manner.  That  this  assumption  was  also  satisfactory 
for  sea  ice  in  the  range  0.  01  <^  «■  <  0.  18  kg/cm2  -  sec  was  demonstrated  by 
Peyton's  data122  which  show  no  pronounced  variation  in  erf  with  changes  in  er 
for  direct  tension  tests  on  sea  ice  from  Point  Barrow,  Alaska,  Peyton122 
then  extended  these  measurements  to  higher  values  of  ir  (~2.  6  kg/cm2  -  sec) 
and  fitted  by  least-squares  a  relation 


.  b  c  c  d  e 
erf  -  a  er  9i  S.  z 


(4.17) 


in  the  linearized  form 


In  er^  =  In  a  +  bln  j  +  c  In  0.  +  d  In  S.  +  e  In  z 


(4.  18) 


to  the  test  results.  Here  a,  ....  e  are  constants,  erf  the  failure  strength, 
er  the  rate  of  stress  application,  0j  the  test  temperature,  Sj  the  ice  salinity, 
and  z  the  position  of  the  sample  in  the  ice  sheet.  The  average  value  of  the 
coefficient  b  was  found  to  be  0.  05  indicating  to  us  that  erf  is  essentially 
independent  of  ir. 

Until  additional  evidence  to  the  contrary  becomes  available,  it  is  also 
probably  satisfactory  to  assume  that  cTf  values  determined  by  ring  tensile 
tests  are  independent  of  ir  in  the  ir  range  usually  encountered  in  these  tests 
(1<  cr  <  3  kg/cm2  -  sec).  An  examination  of  the  results  of  Butkovich34 
supports  this  suggestion,  which  should,  of  course,  be  experimentally  veri¬ 
fied. 

Flexural  strength.  The  effect  of  variations  in  or  on  the  value  of  the 
flexural  strength  as  determined  by  small  beam  tests  has  been  studied  in  con¬ 
siderable  detail  by  Tabata.  **2 » I*e  During  each  constant  deflection-rate 
test  a  plot  of  applied  load  vs  time  was  recorded  and  the  failure  was  classi¬ 
fied  as  elastic  or  "plastic"  depending  upon  whether  or  not  the  plot  was  linear 
or  showed  a  curvature.  The  percentage  of  samples  that  failed  elastically 
increased  from  zero  at  low  values  of  (a-  <  0,2  kg  /cm2  -  sec)  to  100%  at  high 
values  of  cr  (cr^  0.6  kg/cm2  -  sec).  Between  these  two  regions  there  is  a 
mixed  region  in  which  both  elastic  and  plastic  breaks  occur  (Fig.  42).  The 
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boundaries  separating  these  different  "failure  types"  are,  of  course, 
dependent  on  the  test  temperature;  higher  values  of  b  are  required  to 
achieve  elastic  failure  at  temperatures  near  the  melting  point.  It  will, 
in  the  future,  be  interesting  to  separate  the  effects  of  variations  in  sample 
temperature  from  effects  associated  with  variations  in  the  volume  of  brine 
and  air  in  the  ice. 

When  constant  temperature  plots  of  erf  vs  b  are  prepared  and  only 
specimens  that  failed  elastically  are  included,  o-f  appears  to  be  independent 
of  b  in  the  range  0.3  to  8  kg/cm2  -  sec.  If,  however,  both  plastic  and 
elastic  breaks  are  considered  (Fig.  43),  as  cr  increases  there  first  occurs 
a  slight  decrease  (~1  kg/cm2  )  in  and  then  an  increase.  The  minimum 
occurs  in  the  range  of  the  "mixed  region"  between  plastic  and  elastic 
failures.  Tabata*42  suggests  that  the  increase  in  o"£  with  decreasing  or  in 
the  low  b  range  is  due  to  plastic  deformation  of  the  samples  before  breaking. 
The  increase  in  <j-f  with  b  in  the  high  stress  rate  range  is  presumed  to  be 
related  to  the  increase  in  the  static  Young's  modulus  with  increasing  b  at 
b  values  greater  than  3  to  5  kg/cm2  -  sec.  Because  the  samples  used  by 
Tabata  in  these  studies  wer|  obtained  from  very  warm,  new  sea  ice  (20 
to  25  cm  thick)  and  it  was  necessary  to  store  the  ice  for  appreciable  periods 
of  time,  significant  brine  drainage  occurred  prior  to  testing.  How  the 
replacement  ot  the  initial  brine  by  air  affects  the  results  is  currently  un¬ 
known.  It  would  be  desirable  to  repeat  this  careful  series  of  tests  on  sea 
ice  produced  in  the  laboratory  so  that  brine  drainage  prior  to  testing  could 
be  minimized. 

Ta'bata  and  his  co-workers14*’  144,147  > lf*  have  also  studied  the  effect  of 
variations  in  the  manually  applied  average  stress  rate  b  on  the- failure 
strength  as  measured  by  in-situ  cantilever  beam  tests  on  ice  less  than 
30  cm  thick.  Mean  ice  temperatures  were  quite  high,  varying  from  -0.8 
to  -2.  5C.  Ice  salinities  varied  from  4.7  to  13.6%o.  Although  in-situ  tests 
avoid  the  brine  drainage  problem,  the  results  are  difficult  to  analyze  be¬ 
cause  of  the  vertical  variation  in  the  temperature  and  salinity  of  the  ice 
sheet.  The  results  are,  however,  quite  interesting.  Pull-ups,  push-downs 
and  pull-sideways  tests  gave  approximately  the  same  results.  The  value 
of  b  varied  in  the  range  of  0.  1  to  40  kg/cm2  -  sec. 

The  results  of  these  tests  (Fig.  44)  show  that  the  flexural  strength  o-f 
is  independent  of  b  at  values  of  b  .<  1.0  kg /cm*  -  sec.  On  the  other  hand 
when  b  >  1.0  kg/cm2  -  sec,  Of  increases  as  a  linear  function  of  In  b.  That 
the  transition  between  these  two  behavorial  regions  coincides  with  the  change 
from  plastic  to  elastic  deformation  has  been  demonstrated  by  examining 
plots  of  applied  force  vs  deflection  for  individual  beam  tests.  These  plots 
change  from  curved  to  linear  at  values  of  b  >  1  kg/cm2  -  sec.  It  is  reason¬ 

able  to  assume  that  this  transition  to  elastic  behavior  will  occur  at  lower 
values  of  b  for  lower  test  temperatures.  The  slope  also  appears  to  increase 
as  the  ice  temperature  is  lowered. 

Compressive  strength.  The  effect  of  variations  in  b  on  the  value  of  the 
compressive  strength  of  sea  ice  from  Point  Barrow,  Alaska,  has  also  been 
studied  by  Peyton122  • 125  using  eq  4.  17  and  eq  4.  18.  The  range  of  stress 
rates  used  was  0.  06  <  b  <_  8.2  kg/cm2  -  sec.  The  averaue  value  of  b  as 
determined  by  least-squai  cs  was  0.23  (Peyton,  1967).  This  indicates  that 
although  compressive  strength  is  more  dependent  cn  b  tha,t  tensile  strength, 
the  variation  is  still  small.  The  general  increase  in  compressive  strength 
with  an  increase  in  b  is  comparable  to,  although  much  less  prone'  need  than, 
the  increase  in  flexural  strength  with  increasing  j  as  demons. rated  by 
Tabata  and  his  co-workers. 
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Figure  44.  Relation  between  in-situ 
flexural  strength  a-£  and  the  stress  rate 
b  at  -2C  (Tabata,  Fujino  and  Aota). 149 


Figure  45.  Relation  between  compres¬ 
sive  strength  <rf  and  stress  rate  a 
(Peyton). 123 


Peyton123  has  also  presented  some,  results  cf  compressive  strength 
measurements  vs  b  (Fig.  45)  as  determined  on  ice  from  Cook  Inlet,  Alaska. 
These  tests  indicate  a  sharp  decrease  in  &£  with  increasing  b.  It  is  inter¬ 
esting  that  this  decrease  occurs  in  the  same  range  of  b  values  (compare 
Fig.  43,  44  and  45)  in  which  both  Peyton's  Point  Barrow  compressive  data 
and  Tabata's  flexural  data  indicate  an  increase  in  o-f  with  increasing  b. 

Peyton  has  suggested  that  because  the  general  structure  of  the  ice  at  Point 
Barrow  and  Cook  Inlet  is  different,  the  minimum  b  values  of  the  Cook  Inlet 
tests  can  be  interpreted  to  match  the  Barrow  results  at  the  maximum  strength 
value  of  "comparable"  ice.  This  would  then  indicate  a  general  erf  -  f(cr)  re¬ 
lation  with  a  maximum  value  of  q  at  a  ~1  kg/cm2  -  sec.  The  experimental 
evidence  is,  at  present,  not  sufficient  for  definite  conclusions.  An  overall 
theory  of  the  rate  dependence  is  still  required. 

Elastic  modulus.  It  has  already  been  noted  that  the  values  of  E  deter¬ 
mined  by  static  methods  are  consistently  lower  than  when  determined  by 
dynamic  methods.  Considerable  insight  into  this  difference  has  been  pro¬ 
vided  by  the  work  of  Tabata  and  his  associates  on  the  relation  between 
static  E  values  and  b  using  small  sample  beams162  ,  168  and  in-situ  cantilever 
tests. 165’ 166  ’ 167  ’ 169  In  the  beam  tests,  although  more  data  would  be  desirable, 
log  E  appears  to  increase  with  log  cr,  approaching  the  dynamic  value  at  large 
values  of  Jr.162  The  results  of  the  cantilever  tests  are  shown  in  Figure  46. 
Again  a  pronounced  change  in  E  v/ith  <r  was  noted.  The  increase  in  the  value 
of  E  determined  during  the  1965  field  season  over  that  obtained  in  1964  is 
attributed  to  the  increase  in  the  percentage  of  infiltrated  snow  ice  present 
in  1965.  This  interpretation  is  supported  by  the  results  of  dynamic  meas¬ 
urements  (Fig.  39). 


Creep 

Studies  of  the  viscoelastic  properties  of  sea  ice  have  been  performed 
mainly  jy  Tabata155  • 158  >16J  and  Tabata  and  Ono.160  They  have  utilized  both 
small  sample,  and  in-situ,  test  procedures.  As  Tabata  has  pointed  out,158 
the  in-situ  procedures  avoid  brine  drainage  problems  that  would  be  en¬ 
countered  when  small  samples  are  stressed  for  long  periods  of  time. 
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Figure  46.  Relation  between  elastic  modulus 
E  and  stress  rate  <r  at  - 1 .  5C  as  determined  by 
in-situ  beam  te sts  (  Tabata,  Fujino  and  Aota)^89 


However,  as  mentioned  earlier,  in-situ  tests  are  not  without  their  dis¬ 
advantages.  In  particular  in  long-time  in-situ  tests  the  presence  of  warm 
sea  water  in  the  channel  cut  to  produce  the  beam  will  modify  the  vertical 
temperature  profile  in  the  beam.  Because  of  this  problem  and  ice  formation 
in  the  channels,  Tabata  limited  his  in-situ  tests  to  time  periods  of  less  than 
two  hours. 

The  in-situ  tests  were  of  two  types.  In  one  a  standard  cantilever  was 
prepared  and  the  load  applied  at  the  end  of  the  cantilever.  In  the  other  a 
rectangular  beam  was  cut  which  was  still  fixed  to  the  ice  sheet  at  both  ends. 
The  load  was  applied  to  the  center  of  the  beam.  In  both  cases  the  deflections 
of  the  beams  were  measured  relative  to  the  neighboring  ice.  A  typical 
deflection-time  plot  of  the  results  of  an  in-situ  test  are  shown  in  Figure 
47155.158  Upon  application  of  the  load  we  note  an  instantaneous  elastic 
deformation  (AB),  a  transient  creep  portion  (BC)  and  a  linear  steady  state 
creep  portion  (CD).  When  the  load  is  removed  the  beam  experiences  an 
immediate  elastic  contraction  (DE)  and  then  with  the  passage  of  time 
gradually  approaches  the  value  of  the  permanent  deformation  G  (note  that 
EG  ~BB').  The  time-dependent  portions  of  both  the  loading  and  unloading 
curves  can  be  shown  to  be  to  a  good  approximation  exponential  functions  of 
time.1*8  Results  similar  to  the  loading  portion  of  this  curve  (AD)  have  also 
keen  obtained  by  subjecting  small  sea  ice  samples  to  constant  compressive 
stresses  in  the  laboratory.1*8’180 

Such  a  deflection-time  behavior  can  be  analyzed  by  a  rheological  model 
composed  of  a  Maxwell  unit  and  a  Voigt  unit  connected  in  series  (Fig.  47). 
The  elastic  spring  constants  and  the  viscous  d^sh-pot  constants  associated 
with  these  models  are  E[  ,  E2  and  r| x  ,  r|2  respectively.  In  addition  the 
relaxation  time  Ti  and  the  retardation  time  t2  (ti  =  ru  /Et  ,  t2  =  t^2  /E2  ) 
can  be  obtained.  A  summary  of  the  results  of  these  studies  is  given  in 
Table  II.  It  should  be  noted  that  the  values  of  r| 2  ,  r|2  and  r|  in  Table  II  are 
.related  to  the  specific  tests  under  discussion,  not  to  a  pure  shear  test  as 
is  the  true  physical  constant  t|. 
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Figure  47.  Deflection-time  curve  for  sea 
ice. 155,158 


Tile  in-situ  results  show  appreciable  changes  depending  upon  whether 
the  beam  is  deflected  vertically  or  horizontally.  There  currently  are 
insufficient  in-situ  tests  to  ntudy  any  possible  correlations  between  the 
viscoelastic  parameters  and  the  state  of  the  ice.  With  the  exception  of  the 
values  of  n  ,  the  small  sample  results  are  of  the  same  order  of  magnitude 
as  the  values  from  the  in-situ  "ests.  There  also  does  no',  appear  to  bv  any 
pronounced  variation  in  the  small- sample  elastic  or  viscosity  constants 
with  the  state  of  the  sea  ice.  However,  the  values  of  E2  ,  t),  and  ti2  decrease 
as  the  angle  $  between  the  axis;  of  the  cylindrical  test  specimen  and  the  ice 
sheet  surface  changes  from  90  to  0  degrees.  Tabata1'8  has  suggested  that 
the  low  viscosity  when  <f>  =  0  degrees  i?  caused  by  slip  between  the  platelets 
that  compose  the  crystals  of  sea  ice.  It  is  also  possible  that  some  of  the 
difference  between  the  values  of  the  constants  at  <t>  =  0  and  90  degrees  could 
be  produced  by  differences  in  the  volume  of  infiltrated  snow  ice  present  in 
the  sample.  For  many  practical  purposes,  the  utilization  of  the  four  con¬ 
stants  associated  with  the  combined  Maxwell- Voigt  model  is  inconvenient. 
Therefore,  Tabata  has  also  used  a  simple  Maxwell  model  to  study  his  test 
results.  The  values  of  the  Maxwell  model  constants  F.  and  t  are  also  given 
in  Table  II . 

In  a  later  paper,1*3  Tabata  has  analyzed  the  results  of  time  dependent 
deflection  studies  cn  fixed  end  beams  by  using  an  elastic  approach  and  taking 
the  tension  of  the  neutral  surface  into  account.  He  obtained  values  of  E  from 
the  initial  instantaneous  deflection.  E  is  shown  to  decrease  as  a  linear 
function  of  salinity  and  presumably  also  of  brine  volume  as  indicated  by 
dynamic  tests.  As  was  anticipated  in  the  discussion  of  the  relative  magnitudes 
of  static  and  dynamic  values  of  E,  these  E  values  are  considerably  lower  than 
the  dynamic  values.  The.  creep  rate  was  studied  by  replacing  E  in  the  elastic 
deflection  equation  with  a  time  dependent  viscosity  term  r|(t)  and  evaluating 
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Figure  48.  Typical  examples  of  the  rela¬ 
tion  between  the  coefficient  of  viscosity  r) 
(t)  of  sea  ice  and  the  time  of  loading  as 
determined  from  in-  situ  fixed-ended  beams 
(Tabata).163 


Figure  49.  Coefficient  of  viscosity^ 
vs  ice  temperature  for  specimens  of 
different  salinities.  Curves  1  and  2 
based  on  flexural  studies,  curve  3 
based  on  torsion  (Nazintsev). 103 


T|(t)  as  a  function  of  t.  Some  representative  results  are  shown  in  Figure  48 
and  suggest  a  power  law  with  an  average  power  of  0.  67.  This  indicates  that 
the  relation  between  strain  e  and  time  t  is 


where  k  is  a  constant.  Relation  4.  19  is  similar  to  the  empirical  Andrade 
creep  law.  Large  scale  tank  tests64  have  also  produced  similar  deflection 
vs  time  profiles  to  those  observed  by  Tabata. 

Nazintsev103  has  also  used  small- sample  tests  (rods  in  torsion  and  beams 
in  flexure)  to  study  the  time-dependent  deformation  of  sea  ice  in  the  tempera¬ 
ture  range  -4  to  -28C.  At  constant  stress  he  found  that  only  after  3  days  was 
a  consistent  relation  observed  between  strain  rate  £  and  temperature. 

Figure  49  shows  some  of  his  better  results.  In  this  figure  curves  1  and  2  were 
determined  in  flexure  and  curve  3  in  torsion.  A  pronounced  variation  ot  vis¬ 
cosity  with  salinity  and  temperature  is  suggested.  When  the  present  authors 
attempted  to  analyze  his  complete  results  (his  Table  1,  p.  52)  in  terms  of 
brine  volume,  no  consistent  relation  was  noted.  Considering  the  fact  that  the 
majority  of  Nazintscv's  tests  were  performed  on  young  or  annual  sea  ice,  the 
extremely  low  salinities  reported  suggest  that  significant  brine  drainage 
occurred  during  testing.  This  drainage  may  explain  the  lack  of  correlation 
between  t)  and  vjj.  Although  there  was  a  complete  overlap  of  the  values  de¬ 
termined  by  torsion  and  flexure  in  the  low  viscosity  range,  viscosity  values 
higher  than  50  dy-sec/cm2  were  observed  only  in  flexure.  When  Nazintsev 
attempted  to  analyze  his  results  in  terms  of  a  "Glen-type"  power  law  between 
strain  rate  and  shear  stress,  he  observed  that  the  power  n  was  roughly  1.6  as 
compared  with  a  value  of  5.5  for  glacier  ice.  One  should  consider,  however 
that  the  use  of  the  power  law  itself  is  empirical.  The  relation  of  strain  rate 
to  stress  should  be  linear  at  low  stress  levels. 
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Figure  50.  Stress-rupture  behavior  of  nat¬ 
ural  sea  ice  at  -20C  (range  -10C  to  -28C). 
Dashed  lines  indicate  the  95%  confidence 
limits  (Kingery  and  French).77 


The  time-dependent  stress-rupture  behavior  of  sea  ice  has  also  been 
studied  by  Kingery  and  French77  using  small  simply- supported  beams.  The 
results  (Fig.  50)  indicate  a  linear  relation  between  stress  and  log  (rupture 
time).  There  is  an  appreciable  scatter  in  the  results  which  is  undoubtedly  in 
large  part  produced  by  the  wide  range  of  temperatures  (-10  to  -29C)  at  which 
the  tests  were  performed.  It  would  be  interesting  to  study  changes  in  the 
creep-rupture  behavior  of  sea  ice  with  changes  in  the  state  of  the  ice  at  the 
time  of  testing. 


5.  PLATE  CHARACTERISTICS 

There  is  a  wide  range  of  applied  problems,  the  solution  of  which  depends 
upon  a  knowledge  of  the  mechanical  properties  of  sea  ice.  It  is  beyond  the 
scope  of  this  review  to  discuss  these  in  detail.  A  comprehensive  selection  of 
papers  dealing  with  these  problems  is,  however,  included  ir.  the  reference 
list.  For  a  representative  sampling  of  these  papers  see  Hertz  (1884), 

Bernshtein  ( 1 929),  Shul'man  (1946),  Wyman  (1950),  Nevel(1958,  1966) 

Meyerhof  (I960),  Assur  ( 1  961a),  and  Cutcliffe,  Kingery  and  Coble  ( 1 963). 

All  of  these  problems  have  certain  features  in  common:  it  is  usually  necessary 
to  be  able  to  estimate  effective  values  of  flexural  strength  crp  elastic  modulus 
E,  and  Poisson's  ratio  p  for  the  ice  sheet  as  a  whole.  Using  these  values  it 
is  then  necessary  to  calculate  such  composite  properties  as  the  flexural  rigidity, 
the  action  radius,  the  section  modulus,  the  moment  of  inertia,  and  the  position 
of  the  neutral  axis.  In  the  past  the  considerable  vertical  variation  in  the  phys¬ 
ical  properties  of  the  sea  ice  has  usually  been  ignored  and  the  ice  sheet  has 
been  assumed  to  be  homogeneous.  In  concluding  this  paper  we  would  like  to 
show  how  the  vertical  variation  in  the  physical  properties  of  sea  ice  sheets 
can  readily  be  incorporated  into  the  treatment  of  problems  that  deal  with  the 
ice  sheet  as  a  whole.  22 


_  .  -  -  -  -  ■■■—> 


62 


THE  MECHANICAL  PROPERTIES  OF  SEA  ICE 


r 


Figure  51.  Series  of  schematic 
salinity  profiles  for  sea  ice  of 
various  thicknesses. 


Figure  52.  Series  of  schematic  sa¬ 
linity  profiles  for  sea  ice  100,  200, 
and  300  cm  thick.  The  .ow  salinity 
values  at  the  top  of  the  200  cm  pro¬ 
file  indicate  that  this  ice  has  been 
through  a  period  of  summer  melt. 


As  has  been  demonstrated,  the  main  variable  controlling  the  vertical 
variation  in  the  physical  characteristics  of  sea  ice  is  the  vertical  variation 
of  the  volume  of  brine  in  the  ice  vjj.  The  value  of  r‘  any  position  in  the 
ice  sheet  is,  in  turn,  determined  via  the  phase  relations  by  the  ice  tempera¬ 
ture  61  and  salinity  Sj  at  that  point.  Fortunately  ice  salinity  profiles  go 
chrough  a  systematic  series  of  changes  as  ice  growth  progresses.  Repre¬ 
sentative  values  for  salinity  profiles  can  be  obtained  by  examining  the  pub¬ 
lished  literature97’186  * 145  and  are  shown  schematically  in  Figures  51  and 
52.  These  profiles  can  be  adjusted  so  that  the  variation  in  salinity  is  a 
smooth  function  of  the  ice  thickness  h.  If  the  ice  temperature  profile  is 
known,  a  brine  volume  profile  can  be  calculated  by  using  the  brine  volume 
tables  published  by  Assur.17  Recently  Frankenstein  and  Garner51  have  cal¬ 
culated  a  linear  expression  based  upon  these  tables 

.  Si  T  49.185  ,  „  1  lc  ,, 

Vb  HP  —  +  °'532  J 

where  is  expressed  in  %o,  0t  in  degrees  C  and  in  absolute  ratios.  The 
approximate  relation  holds  in  tne  temperature  range  -  0.  5  >  6j  >  -  22.  9C 

with  a  correlation  coefficient  of  0.9995.  One  can  estimate  the  ice  tempera¬ 
tures  from  the  ice  surface  temperature  Qo  by  assuming  a  linear  profile  which 
in  many  cases  is  a  good  approximation.93 

One  must  subtract  the  thickness  of  the  strengthless  skeleton  layer 
(hsk  ~2.  5  cm)  from  the  total  ice  thickness.  This  correction  becomes  quite 
important  when  the  properties  of  thin  ice  sheets  are  considered.  Now  express 
position  in  the  ice  sheet  in  terms  of  a  transformed  dimensionless  distance  'z 
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where  z  varies  from  a  value  of  0  at  the  bridging  layer  at  the  top  of  the 
skeleton  layer  to  a  value  of  1  at  the  upper  surface  of  the  ice  sheet 


z 


h  -  hsk  -  Z 

~mr 


sk 


(5.2) 


Here  h  is  the  total  ice  thickness  including  the  skeleton  layer  and  z  is  the 
vertical  distance  to  any  position  in  the  ice  sheet  measured  from  the  upper 
surface  of  the  sheet.  If  a  linear  approximation  is  used  for  the  temperature 
profile,  the  value  of  0^  at  any  z  is 

0.  =  0k  -  A0  Z  (5.3) 

where 

A0  =  Go  0b  (5.4) 

and 

«b  *  «»  -  «>w  -  (t)  <5-5> 


where  0j-,  is  the  temperature  at  the  bridging  layer  and  0^,  the  water  tempera¬ 
ture.  If  a  linear  temperature  profile  and  a  constant  ice  salinity  are  assumed, 
it  is  possible  over  a  limited  temperature  range  to  express  the  brine-volume 
profile  analytically  in  terms  of  the  salinity  and  density  of  the  brine  and  ice. 
Because  of  these  restrictions,  Assur22  has  also  suggested  the  expression  of 
and  s/vjj  as  polynomials  in'?  and  utilize  the  relations 


(5.6) 


(5.7) 


to  express  the  variation  of  E  and  o^with  v,  ,  We  expand  here  these  ideas 
considering  that  at  high  brine  volumes  <j-£  does  not  approach  the  value  of 
~0.25  suggested  by  the  linear  extrapolation  of  <j-£  determinations  performed 
on  cold  (vvy,  <  0.40)  sea  ice  (see  p.  36).  Instead  <rf  remains  roughly 
constant  at  vvfc  >  0.40  (see  Fig.  26,  31,  32).  As  discussed  earlier,  these 

are  still  just  indications  inasmuch  as  the  problem  of  plastic  stiess  relief  for 
ring  tensile  samples  at  high  temperatures  has  not  been  resolved.  Although 
the  data  available  on  the  behavior  of  E  at  large  brine  volumes  are  very 
limited,  it  is  reasonable  that  E  must  behave  in  a  similar  manner:  if  sea 
ice  possesses  strength  at  values  of  >  0.25,  it  also  must  have  resistance 
to  strain  and  therefore  a  value  of  E.  On  p.  22  and  23  a  method  was  suggested 
for  handling  the  observed  variation  in  g-f  by  using  eq  3.40  and  setting  v0  =  1 
in  that  equation.  This  gives 


TO 


-  ^  (2  -  ^) 


1 


(5.8) 
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where  -2  represents  the  limiting  slope  as  \/vb  approaches  zero.  A  similar 
relation  can  be  constructed  for  E  in  such  a  way  that  it  satisfies  the  available 
data  in  the  low  brine  volume  range 


=  1  - 


H-3vb). 


Note  that  eq  5.8  and  eq  5.9  are  both  zero  when  =  1  and  approach  eq  5.7 
and  eq  5.  6  respectively  with  Vo  =  0.  25  at  small  values  of  v^.  This  can  be 
verified  by  approximating  eq  5.8  and  eq  5.9  with  their  first  derivatives  at 
vb  -0. 

It  is  convenient  to  express  eq  5.  8  and  eq  5.  9  in  terms  of  polynomials 


(V  rv 


=  a0  +  a,  z  +  a2  z 2  +  a3  z  3  +  . . . 


(5.10) 


— —  =  bo  +  bj  z  +  \>i  z  ^  +  b3  z  ^  +  #  # , 

<r0 


(5.11) 


Assur22  shows  how  to  calculate  the  coeffi  cients  simply  upon  the  basis 
of  values  at  selected  levels.  The  average  values  of  E  and  for  the  whole 
ice  sheet  excluding  the  skeleton  layer  now  becomes 

l 

^  ■£-  dz  =  a0+-i^-  +  -^+-^f+...  (5.12) 


i  dz  =  bo  +  h  +  +  £+. 

o'  o  J  <ro  2  3  4 


(5.13) 


Note  that  in  this  approach  it  is  always  possible  to  integrate  over  the  complete 
thickness  (z  =  0  to  z  =  1)  of  the  ice  sheet  i.  e.  the  problem  of  the  possible 
occurrence  of  a  zero  strength  or  E  layer  (v^  >  v0  in  eq  5.6  or  eq  5.  7)  in  the 
ice  sheet  is  avoided. 

The  average  strength  of  the  profile  according  to  eq  5.  13  becomes  im¬ 
portant  in  some  engineering  problems,  such  as  forces  exerted  on  vertical 
piles.  Normally,  however,  one  deals  with  bending  for  which  n.ore  sophis¬ 
ticated  calculations  are  required. 

As  a  result  of  the  vertical  variation  in  E,  the  neutral  axis  is  no  longer 
in  the  center  of  the  ice  sheet.  Finding  its  location  is  similar  to  finding  the 
center  of  gravity  of  a  cross  section  with  a  width  in  each  layer  reduced 
according  to  the  ratio  E/E0.  Defining  the  nit  moment  as 

**n  =  I  "n  TT  d?  {5‘14> 

o 


(5.14) 
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where  the  superscript  ^indicates  that 
the  parameter  is  calculated  using  the^ 
reduced  cross  section  (see  Fig.  53),  zn 
the  non-dimensional  distance  to  the 
neutral  axis  measured  from  the  bridging 
layer  at  the  bottom  of  the  ice  sheet  is 


(5.  15) 


fit!  = 


(5.16) 


Figure  53.  Relative  width  of 
effective  cross-section(E/E0) 
vs  dimensionless  ice  thick- 


and  fit0  from  eq  5.  12.  As  eq  5.15  clearly 
shows,  the  position  of  the  neutral  axis 
will  shift  with  the  coefficients  at  ,  a2  ,  ... 
which  in  turn  are  determined  by  the  salinity 
and  temperature  profiles. 

To  obtain  D,  the  flexural  rigidity  of 
the  ice  sheet,  we  calculate  the  moment 
of  inertia  f  of  our  reduced  cross  section 
using  the  general  relation 


T  =  fit. 


g,* 

‘  tC 


(5.17) 


E0  rhb3 


Eo  h.  > 


(5.18) 


»  %  +  +  if  +  *  +  - 


(5.19) 


The  term 


h.  =  h  -  h  . 
b  sk 


(5.20) 


appears  when  the  transformation  is  made  back  to  dimensional  units  from  the 
non-dimensional  z  scale.  The  action  radius  l  can  easily  be  obtained  from 
eq  5.  18  by  using  eq  1.5.  Thus  the  deflections  and  moments  for  all  pre¬ 
viously  solved  classical  plate  problems  may  be  utilized  provided  t  for  sea 
ice  is  defined  using  eq  5.  18. 

It  now  remains  to  define  the  flexural  strength  for  use  ir  ..Tate  equa¬ 
tions.  The  stress  in  the  bottom  fiber  is 


(5.21) 


Now  <rb  and  Eb  or  and  Ej  can  be  calculated  from  eq  5.8  and  eq  5.9  by 
introducing  vb  for  the  bottom  or  the  top,  respectively.  Note  that  the 
product  <rbE0 /Eb  in  eq  5.25  is  defined  by  the  state  of  the  bottom  fiber 
which  is  relatively  constant.  5?^  on  the  other  hand,  is  highly  weather 
dependent.  In  the  top  fiber  (eq  5.  26)  both  o-tE0  / Et  and  are  weather 
dependent. 

We  started  this  review  with  practical  problems,  discussed  the  chemical, 
crystallographic,  mechanical  and  physical  aspects  involved  and  in  conclusion 
showed  how  to  utilize  this  knowledge  to  solve  practical  problems. 
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ice  strength  (phase  relations,  air  volume,  and  structural  considerations), 

(2)  theoretical  considerations  (strength  models,  air  bubbles  and  salt  rein¬ 
forcement,  and  interrelations  between  growth  conditions  and  strength),  (3) 
experimental  results  (tensile,  flexural,  shear,  and  compressive  strength, 
elastic  modulus,  shear  modulus  and  Poisson's  ratio,  time  dependent  effects, 
and  creep),  and  (4)  plate  characteristics.  The  paper  includes  a  review  of 
problems  in  sea  ice  investigations,  relates  the  chemical,  crystallographic, 
mechanical,  and  physical  aspects  involved,  and  concludes  by  showing  how  to 
utilize  this  knowledge  to  solve  practical  problems. 
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